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Abstract 

We study highly oscillating solutions to a class of weakly well-posed hyperbolic initial boundary 
value problems. Weak well-posedness is associated with an amplification phenomenon of oscillating 
waves on the boundary. In the previous works [CGW14, CW14], we have rigorously justified a weakly 
nonlinear regime for semilinear problems. In that case, the forcing term on the boundary has amplitude 
O(e^) and oscillates at a frequency 0(I/e). The corresponding exact solution, which has been shown to 
exist on a time interval that is independent of e € (0,1], has amplitude 0{e). In this paper, we deal with 
the exact same scaling, namely 0(e^) forcing term on the boundary and 0(e) solution, for quasilinear 
problems. In analogy with [CGM03], this corresponds to a strongly nonlinear regime, and our main 
result proves solvability for the corresponding WKB cascade of equations, which yields existence of 
approximate solutions on a time interval that is independent of £ € (0,1]. Existence of exact solutions 
close to approximate ones is a stability issue which, as shown in [GGM03], highly depends on the 
hyperbolic system and on the boundary conditions; we do not address that question here. 

This work encompasses previous formal expansions in the case of weakly stable shock waves [MR83] 
and two-dimensional compressible vortex sheets [AM87]. In particular, we prove well-posedness for the 
leading amplitude equation (the “Mach stem equation”) of [MR83] and generalize its derivation to a 
large class of hyperbolic boundary value problems and to periodic forcing terms. The latter case is 
solved under a crucial nonresonant assumption and a small divisor condition. 
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1 Introduction 


1.1 General presentation 

This article is devoted to the analysis of high frequency solutions to quasilinear hyperbolic initial 
boundary value problems. Up to now, the rigorous construction of such solutions is known in only a 
few situations and highly depends on the well-posedness properties of the boundary value problem one 
considers. In the case where the so-called uniform Kreiss-Lopatinskii condition is satisfied, the existence 
of exact oscillating solutions on a fixed time interval has been proved by one of the authors in [Wil02] , see 
also [Wil96] for semilinear problems. The asymptotic behavior of exact solutions as the wavelength tends 
to zero is described in [CGWll] for wavetrains and in [CW13] for pulses. The main difference between 
the two problems is that in the wavetrains case, resonances can occur between a combination of three 
phases, giving rise to integro-differential terms in the equation that governs the leading amplitude of the 
solution^. Resonances do not occur at the leading order^ for pulses, which makes the leading amplitude 
equation easier to deal with in that case. 

In this article, we pursue our study of weakly well-posed problems and consider situations where 
the uniform Kreiss-Lopatinskii condition breaks down. Let us recall that in that case, high frequency 
oscillations can be amplified when reflected on the boundary. As far as we know, this phenomenon was 
first identified by Majda and his collaborators, see for instance [MR83, AM87, MA88] in connection with 
the formation of specific wave patterns in compressible fluid dynamics. Asymptotic expansions in the 
spirit of [AM87] are also performed in the recent work [WY14]. In various situations (depending on the 
scaling of the source terms and on the number of phases), these authors managed to derive an equation 
that governs the leading amplitude of the solution. Solving the leading amplitude equation in [MR83] and 
constructing exact and/or approximate oscillating solutions was left open. As far as we know, the rigorous 
justification of such expansions has not been considered in the literature so far. The present article follows 
previous works where we have given a rigorous justification of the amplification phenomenon: first for 
linear problems in [CGIO], and then for semilinear problems in [CGW14, GW14]. These previous works 
considered either linear problems, or a weakly nonlinear regime of oscillating solutions for which the 
existence and asymptotic behavior of exact oscillating solutions can be studied on a fixed time interval. 

The regime considered in [MR83, AM87], and that we shall also consider in this article, goes beyond 
the one considered in [GGW14, CW14]. In analogy with [GGM03], this regime will be referred to as that 
of strong oscillations. We extend the analysis of [MR83, AM87] to a general framework, not restricted to 
the system of gas dynamics, and explain why the problem of vortex sheets considered in [AM87] and the 
analogous one in [WY14] yield a much simpler equation than the problem of shock waves in [MR83]. We 
also clarify the causality arguments used in [MR83, AM87] to discard some of the terms in the (formal) 
asymptotic expansion of the highly oscillating solution. We need however to make a crucial assumption in 
order to analyze this asymptotic expansion, namely we need to assume that no resonance occurs between 
the phases. This is no major concern for pulses because interactions are not visible at the leading order, 
and this may be the reason why this aspect was not mentioned in [MR83]. Resonances can have far worse 
consequences when dealing with wavetrains, and what saves the day in [AM87, WY14] is that there are 
too few phases to allow for resonances. This explains why the amplitude equation in [AM87] and the 
corresponding one in [WY14] reduce to the standard Burgers equation. When the system admits at least 

^This is not specific to the boundary conditions and is also true for the Cauchy problem in the whole space, see, e.e., 
[HMR86, JMR95]. 

^Interactions between pulses associated with different phases need to be considered only when dealing with the construction 
of correctors to the leading amplitude. 
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three phases (two incoming and one outgoing), and even in the absence of resonances, the amplification 
phenomenon gives rise, as in [MR83], to integro-differential terms in the equation for the amplitude that 
determines the trace of the leading profile. We refer to the latter equation as “the Mach stem equation”, 
and show how it arises more generally in weakly stable (WR class^) hyperbolic boundary problems with 
a strongly nonlinear scaling, both in the wavetrain setting, where the equation we derive appears to be 
completely new, and in the pulse setting, where the equation coincides with the one derived in [MR83]. 

Our main results establish the well-posedness of the Mach stem equation in both settings, and then 
use those solvability results to construct approximate highly oscillating solutions on a fixed time interval 
to the underlying hyperbolic boundary value problems. In the wavetrain case we are able to construct 
approximate solutions of arbitrarily high order under a crucial nonresonant condition and a small divisor 
condition; in the pulse setting we construct approximate solutions up to the point at which further 
“correctors” are too large to be regarded as correctors. 

In Appendix B we compute the formal large period limit of the Mach stem equation for wavetrains, and 
find a surprising discrepancy (described further below) between that limit and the Mach stem equation 
for pulses derived in [MR83]. 

The construction of exact oscillating solutions close to approximate ones is a stability issue that is far 
from obvious in such a strong scaling. We refer to [CGM03] for indications on possible instability issues 
and postpone the stability problem in our context to a future work. In any case, it is very likely that no 
general answer can be given and that stability vs instability of the family of approximate solutions will 
depend on the system and/or on the boundary conditions, see for instance [CGM04] for further results in 
this direction. 

The precise Mach stem formation mechanism described in [MR83] for reacting shock fronts comes from 
wave breaking (blow-up) in the ’’Mach stem equation”. The numerical simulations in [MR84] suggest 
that the latter equation displays a similar blow-up phenomenon as the corresponding one for the Burgers 
equation. Our results show that, in a precise functional setting, the Mach stem equation is a semilinear 
perturbation of the Burgers equation, which might suggest that the hint in [MR84] is true. However, the 
semilinear perturbation in the Mach stem equation takes the form of a bilinear Fourier multiplier which 
makes the rigorous justification of a blow-up result difficult. We also postpone this rigorous justification 
to a future work. 

Notation 

Throughout this article, we let A4„^Ar(]K) denote the set of n x N matrices with entries in IK = M or C, 
and we use the notation A4Ar(]K) when n = N. We let I denote the identity matrix, without mentioning 
the dimension. The norm of a (column) vector X € is \X\ := (X*X)^/^, where the row vector X* 
denotes the conjugate transpose of X. If X,Y are two vectors in C^, we let X - Y denote the quantity 
'^jXjYj, which coincides with the usual scalar product in when X and Y are real. We often use 
Einstein’s summation convention in order to make some expressions easier to read. 

The letter C always denotes a positive constant that may vary from line to line or within the same 
line. Dependance of the constant C on various parameters is made precise throughout the text. The sign 
< means < up to a multiplicative constant. 

^The WR class is described in Assumption 1.6 below. 
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1.2 The equations and main assumptions 


In the space domain := {x = (?/, G ^ x M : > 0}, we consider the quasilinear evolution 

problem with oscillating source term: 


( 1 . 1 ) 


dtUe + Y7j=l ('«£;) djUe = 0 , 
Kue\xd=o) = e^G y, 

'^61 G'|t<o ~ 0) 


t < T, X ^ , 

t <T, y £ , 


where the Aj’s belong to AtAr(M) and depend in a C°° way on u in a neighborhood of 0 in R^, 6 is a 
C°° mapping from a neighborhood of 0 in R'^ to R^ (the integer p is made precise below), and the source 
term G is valued in RP. It is also assumed that 6(0) = 0, so that the solution starts from the rest state 
0 in negative times and is ignited by the small oscillating source term G on the boundary in positive 
times. The two main underlying questions of nonlinear geometric optics are: 


1. Proving existence of solutions to (1.1) on a fixed time interval (the time T > 0 should be independent 
of the wavelength e G (0,1]). 

2. Studying the asymptotic behavior of the sequence as e tends to zero. If we let denote an 
approximate solution on [0, T'], T' < T, constructed by the methods of nonlinear geometric optics 
(that is, solving eikonal equations for phases and suitable transport equations for profiles), how well 
does approximate Ue for e small ? For example, is it true that 


lim \\ue - ^ 0? 


The above questions are dealt with in a different way according to the functional setting chosen for 
the source term G in (1.1). More precisely, we distinguish between: 

• Wavetrains, for which G is a function defined on (—oo, Tq] x R'^“^ x R that is 0-periodic with respect 
to its last argument (denoted 6q later on). 

• Pulses, for which G is a function defined on (—oo, Tq] x R'^"^ x R that has at least polynomial decay 
at infinity with respect to its last argument. 


The answer to the above two questions highly depends on the well-posedness of the linearized system 
at the origin: 


( 1 . 2 ) 


dtv + ^i(O) djV = f , t < T, X G R^ , 

< d6(0) • = g, t <T, y £ R'^"^, 

^v,f,g\t<o = 0. 


The hrst main assumption for the linearized problem (1.2) deals with hyperbolicity. 

Assumption 1.1 (Hyperbolicity with constant multiplicity). There exist an integer q > 1, some real 
functions Ai,...,Aq that are analytic on R*^ \ {0} and homogeneous of degree 1, and there exist some 
positive integers vi,... ,iyq such that: 


VC = (6,...,ed)eM''\{0}, 


det 




j -^j ( 0 ) 


(''■ + ^kiO) 




i=i 


k=l 
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Moreover the eigenvalues Ai(^),..., Ag(^) are semi-simple (their algebraic multiplicity equals their geo¬ 
metric multiplicity) and satisfy Ai(^) < • • • < Xq{C) for all ^ \ {0}. 

For reasons that will be fully explained in Sections 2 and 4, we make a technical complementary assump¬ 
tion. 


Assumption 1.2 (Strict hyperbolicity or conservative structure). In Assumption 1.1, either all in¬ 
tegers equal 1 (which means that the operator dt -\- Aj{0) dj is strictly hyperbolic), or 

Ai{u),... ,A(i{u) are Jacobian matrices of some flux functions fi, ■ ■ ■, fd that depend in a way on 
u in a neighborhood of 0 in . In the latter case. Assumption 1.1 holds for all u close to the origin, 
namely for an open neighborhood O o/O G there exist an integer q > I, some real functions Xi,... ,Xq 
that are C°° on Ox M'^\ {0} and homogeneous of degree 1 and analytic in and there exist some positive 
integers ui,... ,Uq such that: 


det 


d 

tI + 


i=i 


— {r + Afc(tt, .^)) , 


k=l 


for u G O, r G M and ^ = {fi,. .. ,f,d) £ \ {0}. Moreover the eigenvalues Xi{u,f),...,Xq{u,f) are 

semi-simple and satisfy Ai(m,^) < • • • < Xq{u,f,) for all u € O, = (^i,... G \ {0}. 

For simplicity, we restrict our analysis to noncharacteristic boundaries and therefore make the following: 

Assumption 1.3 (Noncharacteristic boundary). The matrix A(i(0) is invertible and the Jacobian matrix 
B := d6(0) has maximal rank, its rank p being equal to the number of positive eigenvalues of Arf(O) 
(counted with their multiplicity). Moreover, the integer p satisfies 1 < p < N — 1. 

Energy estimates for solutions to (1.2) are based on the normal mode analysis, see, e.g., [BGS07, 
chapter 4]. We let t — i^ G C and p G denote the dual variables of t and y in the Laplace and 

Fourier transform, and we introduce the symbol 


A(C) :=-i-4rf(0) ^ ^(r-i 7 )/-F^ 7 /jAj( 0 )j , C := (^ - * 7 , € C x ^ 

For future use, we also define the following sets of frequencies: 

S:={(r-Z 7 , 7 ?) gCxM''-1\(0,0) : 7 > o} , S := {c G H : + 7^ + = i| , 

Ho := |(r,7) G M x \ (0,0)| = H n {7 = 0} , Sq := S D Hq. 

Two key objects in our analysis are the hyperbolic region and the glancing set that are defined as follows. 

Definition 1.4. • The hyperbolic region % is the set of all (T,rf) G Hq such that the matrix A(r, 7 ) 

is diagonalizable with purely imaginary eigenvalues. 

• Let G denote the set of all (r, ,^) G M x such that ^ 7^ 0 and there exists an integer /c G {1,... , 7} 
satisfying 

T + Afc(^) = = 0 • 

//vr(G) denotes the projection of G on the first d coordinates (that is, '/r(r, ^) := (t, ^ 1 ,..., ^rf_i) 
for all {t,^)), the glancing set Q is Q := 7 r(G) C Hq. 
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We recall the following result that is due to Kreiss [Kre70] in the strictly hyperbolic case (when all integers 
Vj in Assumption 1.1 equal 1) and to Metivier [MetOO] in our more general framework. 

Theorem 1.5 ([Kre70, MetOO]). Let Assumptions 1.1 and 1.3 he satisfied. Then for all C £ 'B\ Hq, the 
matrix M(C) has no purely imaginary eigenvalue and its stable subspace ]E^(C) ho,s dimension p. Further¬ 
more, E® defines an analytic vector bundle over 'F,\'F,q that can be extended as a continuous vector bundle 
over H. 

For all (r, ry) € Sq, we let E^(r, 77) denote the continuous extension of E® to the point {T,r]). Away from 
the glancing set Q C Eq, E®(C) depends analytically on C, see [MetOO]. In particular, it follows from the 
analysis in [MetOO], see similar arguments in [BGRSZ02, Coull], that the hyperbolic region LL does not 
contain any glancing point, and E^(^) depends analytically on ^ in the neighborhood of any point of Li. 
We now make our weak stability condition precise (recall the notation B := d6(0)). 

Assumption 1.6 (Weak Kreiss-Lopatinskii condition). • For all Q £'F. \ Hq, KerR n E^((^) = {0}. 

• The set T := {(^ € Sq : KeiB nE^(^) yZ {0}} is nonempty and included in the hyperbolic region Li. 

• There exists a neighborhood V of T in T,, a real valued C°° function a defined on V, a basis 

..., Ep{() of E^(C) that is of class C°° with respect to ( £ V, and a matrix P{C) £ GLp(C) 
that is of class C°° with respect to Q £V, such that 

VCgV, B{Ei{C) ... £^^(0) =P(C)diag (7 + iu(C),l,...,l). 

As explained in [GGIO, CGW14, CW14], Assumption 1.6 is a more convenient description of the so- 
called WR class of [BGRSZ02] . Let us recall that this class consists of hyperbolic boundary value problems 
for which the uniform Kreiss-Lopatinskii condition breaks down ”at first order” in the hyperbolic region^. 
This class is meaningful for nonlinear problems because it is stable by perturbations of the matrices Aj{0) 
and of the boundary conditions B. 

Our final assumption deals with the phase ipo occuring in (1.1). 

Assumption 1.7 (Critical phase). The phase ipo in (1.1) is defined by 

To{t,y) ■=zt + v-y, 

with (t, ry) £T. In particular, there holds {T,r]) £ Li. 

The shock waves problem considered in [MR83] enters the framework defined by Assumptions 1.1, 
1.2, 1.3, 1.6 and 1.7 with the additional difficulty that the space domain has a free boundary. The 
vortex sheets problem considered in [AM87] and the analogous one in [WY14] violate Assumption 1.3 
but these problems share all main features which we consider here. We restrict our analysis to fixed 
noncharacteristic boundaries mostly for convenience and simplicity of notation. 

Our main results deal with the existence of approximate solutions to (1.1). This is the reason why 
we only make assumptions on the linearized problem at the origin (1.2), and not on the full nonlinear 
problem (1.1). 

'^Let us also recall that the uniform Kreiss-Lopatinskii condition is satisfied when KerB = {0} for all C G 2, and 

not only for € H \ Ho. 


7 



1.3 Main result for wavetrains 


In Part I, we consider the nonlinear problem (1.1) with a source term G that is ©-periodic with respect 
to its last argument Oq. As evidenced in several previous works, the asymptotic behavior of the solution 
Ug to (1.1) is described in terms of the characteristic phases whose trace on the boundary equals ipQ. We 
thus consider the pairwise distinct roots (and all the roots are real) ,..., to the dispersion relation 


det 


d-l 

i=i 


Aj(0) -|- (jj A,i(0) 


= 0 . 


To each there corresponds a unique integer km € {1,... , such that r -|- — 0. We can 

then define the following real phases and their associated group velocity: 


(1.3) Vm = l,...,M, ipm{t,x) := +i^mXd, := VAfc^(^,w^). 

We let $ := (yji,..., ^pm) denote the collection of phases. Each group velocity is either incoming or 
outgoing with respect to the space domain the last coordinate of is nonzero. This property holds 
because (r, rf) does not belong to the glancing set Q. 

Definition 1.8. The phase tpm is incoming if the group velocity is incoming (d^^Xk^{ri,Um) > 0), 
and it is outgoing if the group velocity Vm is outgoing (d(^^Xk^{ri,uim) < Oj. 

In all what follows, we let X denote the set of indices m € {1,... ,M} such that p>m is incoming, and O 
denote the set of indices m G {1,... ,M} such that ipm is outgoing. Under Assumption 1.3, both I and 
O are nonempty, as follows from [CGIO, Lemma 3.1] which we recall later on. 

The proof of our main result for wavetrains, that is Theorem 1.10 below, heavily relies on the nonres¬ 
onance assumption below. For later use, we introduce the following notation: if 0 < A; < M, we let 
denote the subset of all a G such that at most k coordinates of a are nonzero. For instance Z^d jg 
the union of the sets Z e^, m = 1,..., M, where (ei, ..., bm) denotes the canonical basis of M^. We also 
introduce the notation: 

d d 

(1.4) L(r ,0 :=t/ + ^CjAj( 0 ), L{d) := dt+ ■ 

i=i j=i 

The nonresonance assumption reads as follows. 

Assumption 1.9 (Nonresonance and small divisor condition). The phases are nonresonant, that is for 
all a G Z^ \Z^d^ there holds detL(d(a • <!>)) 7 ^ 0, where a • <I> := amTm- 

Furthermore, there exists a constant c > 0 and a real number 7 such that for all a G Z^ \ Z^d fjiQf 
satisfies am = 0 for all m € O, there holds 


det L(d(Q: • $))| > c |a| . 


Let us note that the small divisor condition is only required for a with nonzero components am which 
correspond to incoming phases. If there is only one incoming phase, then there is no such a with at 
least two nonzero coordinates, and we do not need any small divisor condition. The reason for this 
simplification will be explained in Sections 2 and 3. Our main result reads as follows. 




Theorem 1.10. Let Assumptions 1.1, 1.2, 1.3, 1.6, 1.7, 1.9 he satisfied, let Tg > 0 and consider G G 
C°°((—oo, To]t; X (M/(0Z))0q)) that vanishes for t < 0. Then there exists 0 < T < Tg and 
there exists a unique sequence of profiles fUn )n>o in C°°{{—oo,T]-, x (M/(0Z))^)) that satisfies 

the WKB cascade (4.1), (4.3) below, and h(n\t<o = 0 for all n G N. In particular, each profile Uk has its 
6-spectrum included in the set 


Zf ;= {a G l^m e O, am = O} , 

which means that no outgoing signal is generated at any order. 

Furthermore, if for all integers A^i,A ^2 > 0, we define the approximate solution 

Ni-eN2 

n=0 



then 


d^^.pp,N,,N2 ^ g.^app,iVi,iV2 ^ o^giVi+1) ^ 

i=i 

^ ,2 ^ + 0(e^i+2) , 

^app,iV„iV,|^<0^0^ 


t < T, X £ , 

t <T, y £ , 


where the 0(e'^i'*'^) in the interior equation and in the boundary conditions are measured respec¬ 

tively in theC{{-oo,T\;H^^{W\_))nL°^{{-oo,T\xW\.) andC{{-oo,T\-H^^{Mfi-^))nL°^{{-oo,T\xMfi-^) 
norms. 


Of course, the approximate solutions provided by Theorem 1.10 become interesting only for A^i > 1, 
that is, when the remainder on the boundary becomes smaller than the source term G. 

The spectrum property in Theorem 1.10 is a rigorous justification of the causality arguments used 
in [AM87, WY14]. Theorem 1.10 will be proved in Part I of this article. In Section 2, we shall derive 
the so-called leading amplitude (Mach stem) equation from which the leading profile Hq is constructed. 
Section 3 is devoted to proving well-posedness for this evolution equation. As far as we know, the bilinear 
Fourier multiplier that we shall encounter in this equation had not appeared earlier in the geometric 
optics context and our main task is to prove a tame boundedness estimate for this multiplier. Section 4 
is devoted to the construction of the correctors lAn, n > 1, and to completing the proof of Theorem 1.10. 
We refer to Appendix A for a discussion of the two-dimensional isentropic Euler equations, with specific 
emphasis on Assumption 1.9. 


1.4 Main result for pulses 

We keep the same notation (1.3) for the phases, but now consider the nonlinear problem (1.1) with a 
source term G that has ’’polynomial decay” with respect to its last argument 6q. This behavior is made 
precise by introducing the following weighted Sobolev spaces: 

:= |u G x M^) : 0^ d^gU £ L‘^{R'^) if a + |/3| < fcj . 

Our second main result reads as follows. 
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Theorem 1.11. Let Assumptions 1.1, 1.2, 1.3, 1.6, 1.7 be satisfied. Let denote the smallest integer 
satisfying kQ > {d+ l)/2, and let Kq, Ki be two integers such that Kq > 8 + {d + 2)/2, Ki — Kq >2ko + 2. 
Let Tq > 0 and consider 

that vanishes for t < 0. Then there exists 0 < T < Tq and there exist profiles Uq,Ui,IA 2 vanishing in t < 0 
that satisfy the WKB cascade (5.2), (5.4) below. If we define the approximate solution 


u. 


app 


{t,x) := U,x, 


n=0 


Toit,y) xd 


then 


dtuT^ + Y^AjiuT^) = 0(8^), 

J=1 


t < T, X € 


b{uT^\x,i=o) = e'^G {t,y, 

ur%<o = 0, 


To{t,y) 


+ 0(e^), t<T,y€ 




pd—1 


where the O(e^) in the interior equation and in the boundary conditions are measured respectively in the 
C((-oo,r];L2(M^)) nL°°((-oo,T] x M^) and C{{-oo,T]- L'^{R<^-^)) n L°°{{-oo,T] norms. 

The exact regularity and decay properties of the profiles are given in (5.9). 


The approximate solution provided by Theorem 1.11 is constructed, as in [MR83], according to the 
following line of thought: we expect that the exact solution to (1.1) admits an asymptotic expansion 
of the form 


Ue ^ s'^e'^Uki t,x, 
k>0 


To{t,y) Xd 
e e 


that is either finite up to an order K >2, or infinite. We plug this ansatz and try to identify each profile 
Uk. The corrector eUi is expected to be negligible with respect to Uq, and so on for higher indices. Hence 
the identification of the profiles is based on some boundedness assumption for the correctors to the leading 
profile. Of course, such assumptions have to be verified a posteriori when constructing Uq, Ui and so on. 
For instance, Theorem 1.10 is based on the assumption that one can decompose with profiles in 
and we give a rigorous construction of such profiles for which the corrector 




t,X, 


<I>(t, x) 


is indeed an 0(e^+”') in L°°{{—oo,T] x M^^). 

In Sections 5 and 6, we give a rigorous construction of the leading profile Kq and of the first two 
correctors Ui, IA 2 that satisfy all the boundedness and integrability properties on which the derivation 
of the leading amplitude equation relies. In particular in section 5 we explain why, assuming that the 
hrst and second correctors Ui, U 2 satisfy some boundedness and integrability properties in the stretched 
variables {t,x,6o,(,d), the leading profile Uq is necessarily determined by an amplitude equation that is 
entirely similar to the one in [MR83]. The analysis of Section 5 clarifies some of the causality arguments 
used in [MR83]. This makes the arguments of [MR83] consistent, and one of our achievements is to prove 
in section 6 local well-posedness for the leading amplitude equation derived in [MR83] (which we call the 
Mach stem equation). 
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However, the drawback of this approach is that, surprisingly, it is not consistent with the formal 
large period limit for wavetrains. More precisely, it seems rather reasonable to expect that the pulse 
problem is obtained by considering the analogous problem for wavetrains with a period 0 and by letting 
0 tend to inhnity. In particular, the reader can check that the leading amplitude equations derived in 
[CW13], resp. [CW14], for quasilinear uniformly stable pulse problems, resp. semilinear weakly stable 
pulse problems, coincide with the large period limit of the analogous equations obtained in [CGWll], 
resp. [CGW14], for wavetrains, even though the latter equations include interaction integrals to account 
for resonances. One could therefore adopt a different point of view and first derive the profile equations for 
pulses by considering the limit 0 —)• +oo for wavetrains, and then study the property of the corresponding 
approximate solution. Surprisingly, the two approaches do not give the same leading profile Uq, as we 
shall explain in Appendix B. It seems very difficult at this stage to decide which of the two approximate 
solutions should be the most “physically relevant” since we do not have a nonlinear stability result that 
would claim that the exact solution Ug to (1.1) is close to one of these two approximate solutions on a 
fixed time interval independent of e small enough. The clarification of this surprising phenomenon is left 
to a future work. 
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Part I 


Highly oscillating wavetrains 

2 Construction of approximate solutions: the leading amplitude 

2.1 Some decompositions and notation 

We recall here some useful results from [CGIO] and introduce some notation. Recall that the matrix 
A{t, rf) is diagonalizable with eigenvalues i m = 1,... , M. The eigenspace of A{t, rf) for i coincides 
with the kernel of L{d^pm). 

Lemma 2.1 ([CGIO]). The (extended) stable subspace E®(r, ry) admits the decomposition 

(2.1) E®(t,;^) = ©mexKer L{dLpm ), 

and each vector space in the decomposition (2.1) is of real type (that is, it admits a basis of real vectors). 
Lemma 2.2 ([CGIO]). The following decompositions hold 

(2.2) Ker L{dg^m) = ©™=i ^d(O) Ker , 

and each vector space in the decompositions (2.2) is of real type. 

We let Pi,..., Pm , resp. Qi,... ,Qm, denote the projectors associated with the first, resp. second, 
decomposition in (2.2). Then for all m = 1,..., M, there holds Im L(d(/?m) = Ker Qm- 

Using Lemma 2.2, we may introduce the partial inverse Rm of L(d(pm), which is uniquely determined by 
the relations 

V — 1, . . . , M , Rm 1 Pm j T{^Tm') Rm — I Qm j Pm Rm — 0 , Rm Qm — 0 . 

When the system is strictly hyperbolic, which is the case considered in Sections 2, 3 and most of 
Section 4, each vector space Ker L{d(pm) is one-dimensional and M = N. The case of conservative 
hyperbolic systems with constant multiplicity will be dealt with in Paragraph 4.4. In the case of a strictly 
hyperbolic system, we choose, for all m = 1,..., A^, a real vector that spans Ker L{dipm)- We also 
choose real row vectors ii,... ,i]y, that satisfy 

Vm = l,...,7V, emL{dipm)=0, 

together with the normalization imAd{0)r.m' = ^mm'- With this choice, the partial inverse Rm and the 
projectors Pm, Qm are given by 

yX€C^, RmX= V Trn', PmX = {^mAd{d)X)rm, Qm X = {Im X) Tm ■ 

According to Assumption 1.6, KerR nE^(r, 77 ) is one-dimensional and is therefore spanned by some 
vector e = ^ Span (here we have used Lemma 2.1). The vector space ilE®(r, ry) is 

[p — l)-dimensional and is of real type. We can therefore write it as the kernel of a real linear form 

(2.3) = {X eCP ,bX = d] , 
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for a suitable nonzero row vector b. Eventually, we can introduce the partial inverse of the restriction 
of B to the vector space (r, r/). More precisely, we choose a supplementary vector space of Span e in 
E^(r,^): 

(2.4) E^(r, 77 ) = (r, 77 ) © Span e . 

The matrix B then induces an isomorphism from E^(t, 77) to the hyperplane i?E®(r, 77). 

2.2 Strictly hyperbolic systems of three equations 

For simplicity of notation, we first explain the derivation of the leading amplitude equation in the case of 
a 3 X 3 strictly hyperbolic system. We keep the notation introduced in the previous paragraph, and we 
make the following assumption. 

Assumption 2.3. The phases are incoming and (p 2 is outgoing. 

Assumption 2.3 corresponds to the case p = 2 in Assumption 1.3 (up to reordering the phases). The only 
other possibility that is compatible with Assumption 1.3 is p = 1, and two phases are outgoing. This 
case would yield the standard Burgers equation for determining the leading amplitude (see Paragraph 2.4 
below for a detailed discussion), so we focus on p = 2 which incorporates the main new difficulty. 

Let us now derive the WKB cascade for highly oscillating solutions to (1.1). The solution to (1.1) 
is assumed to have an asymptotic expansion of the form 

(2.5) Ue e Uk (t, x, 

k>o ^ ^ 

where we recall that denotes the collection of phases (pi, (p 2 j y^s), and the profiles Uk are assumed to 
be ©-periodic with respect to each of their last three arguments 6 * 1 , 02 ) ^ 3 - Plugging the ansatz (2.5) in 
(1.1) and identifying powers of e, we obtain the following first three relations for the U^’s (see Section 4 
for the complete set of relations up to any order): 

(a) C{de)Uo = 0, 

( 2 . 6 ) (b) £{de)Ui + L{d)Uo + M{Uo,Uo) = 0, 

(c) Cide)U 2 + Lid)Ui+Mil{o,ldi)+MiUi,Uo) +A'i(Wo,Wo) + A/' 2 (L^o,^^o,Wo) = 0 , 

where the differential operators T, Af,A/i,A /2 are dehned by: 

£{80) := L{d(fm) de^ , 

M{v,w) := djipm (dAj(O) • v) dg^w , 

Ai( 77 , 77 ;) := (dAj(O) • 7 ;) djW , 

J\f 2 {v,v,w) := (d^Aj(O) • {v,v))dg^w. 

The equations (2.6) in the domain (—cx),T] x x (M/0Z)^ are supplemented with the boundary con¬ 
ditions obtained by plugging (2.5) in the boundary conditions of (1.1), which yields (recall B = d6(0)): 

BUq = 0, 

BUi + ^ dH{0) • (UoM) = G{t, y, 0o) , 


( 2 . 8 ) 


(b) 
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where functions on the left hand side of (2.8) are evaluated at = 0 and 6 i = 62 = 03 = Oq- In order to 
get Ue\t<o = 0, as required in (1.1), we also look for profiles Uk that vanish for t < 0. 

The derivation of the leading amplitude equation is split in several steps, which we decompose below 
in order to highlight the (slight) differences in Paragraph 4.4 for the case of systems with constant 
multiplicity. 


Step 1: Uq has mean zero. 

According to Assumption 1.9, the phases (pm are nonresonant. Equation (2.6) (a) thus yields the 
polarization condition for the leading amplitude Uq. More precisely, we expand the amplitude Uq in 
Fourier series with respect to the Om's, and (2.6) (a) shows that only the characteristic modes occur 
in Uq. More precisely, we can write 

3 

(2.9) Uo{t,X,dl, 02 , 63 ) = U.o{t,x) + ^ (Jm{t,X,0m)rm , 

m=l 

with unknown scalar functions am depending on a single periodic variable 9m and of mean zero with 
respect to this variable. 

Let us now consider Equation (2.6)(b), and integrate with respect to 01 , 02 , 03 - Using the expression 

(2.9) of Uq, we obtain 

(2.10) L{d)UQ = 0, 

because the quadratic term M.{Uq,Uq) has zero mean with respect to {0i, 02 , 03 )- Indeed, M{Uq,Uq) splits 
as the sum of terms that have one of the following three forms : 

fmit, ^) , *) f’m , *) m2 ^m\m2 (mi / m2) , 

where fm, rmim 2 are constant vectors (whose precise expression is useless), and each of these terms has zero 
mean with respect to { 01 , 02 , 03 )- Equation (2.10) is supplemented by the boundary condition obtained 
by integrating (2.8) (a), that is, 

( 2 . 11 ) BUq\,^=q = 0. 

By the well-posedness result of [Cou05], we get Uq = 0. The goal is now to identify the amplitudes Um’s 
in (2.9). 


Step 2: Uq has no outgoing mode. 

We first start by showing a 2 = 0. We first integrate (2.6) (b) with respect to { 01 , 03 ) and apply the 
row vector £2 (which amounts to applying Q 2 ), obtaining 


i2L{d){a2r2) + h j~j~ M{Uo,Uq) d0id0:^ =0. 


Since there is no resonance among the phases, integration of the quadratic term A4{Uq,Uq) with respect 
to { 01 , 03 ) only leaves the self-interaction term a 2 902^^2, and the classical Lax lemma [Lax57] for the linear 
part® £2 L{d){a 2 r 2 ) gives the scalar equation 


dta2 + V2 • Va;(T2 + C2 Cr2 5^2Cr2 = 0 , 


C2 := 


djip 2 h (dAj(O) -r 2 )r 2 
^2 r2 


®In fact, £2 L{d){- r 2 ) equals £2 r 2 times the transport operator dt +^2 ■ Va,, and £2 r 2 does not vanish. 
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Since V 2 is outgoing and (T 2 vanishes for t < 0, we obtain (T 2 = 0. 

The above derivation of the interior equation for cj 2 can be performed word for word for the other 
scalar amplitudes cti,c 73 , because A4{Uq,IJq) also has zero mean with respect to { 01 , 62 ) and { 62 , 6 ^). We 
thus get 


( 2 . 12 ) 


T ^x^m T Cm dm 0 , TTl 1)3, Cm 


dj^m (cl^j(O) • Cm) Cj) 


but we now need to determine the trace of dm on the boundary {x^ = 0}. 

Since only di, ds appear in the decomposition (2.9), the leading amplitude Uq takes values in the stable 
subspace IE^(r, ?]) (Lemma 2.1), and the boundary condition (2.8) (a) yields 


di(t, y, 0, Oq) ri = a{t, y, Oq) ci , d 3 (t, y, 0, Oq) = a{t, y, Oq) 63 , 


for a single unknown scalar function a of zero mean with respect to its last argument Oq (recall that 
e = Cl + 63 spans the vector space ]E^(r, y) D Ker B). 


Step 3: Ui has no outgoing mode. 

The derivation of the equation that governs the evolution of a comes from analyzing the equations for 
the first corrector Ui. Since (2.6)(c) is more intricate than the corresponding equation in [CGW14], the 
analysis is starting here to differ from what we did in our previous work [CGW14]. The first corrector Ui 
reads 

3 

Ul{t,X, 01,62,03) = U^{t,x) + ^ U^{t,X, 0 m) +Wi“(t,X,6'i,6»2,6'3) , 

m=l 

where Ui represents the mean value with respect to { 61 , 02 , 63 ), each ZY™ incorporates the 0m-oscillations 
and has mean zero, and the spectrum of the ’’noncharacteristic” part is a subset of \ Z^h f;Q 
the nonresonant Assumption 1.9. More precisely, is obtained by expanding (2.6) (b) in Fourier series 
and retaining only the noncharacteristic modes Z^ \ Prom the expression (2.9) of Uq (recall ^ 0 
and d 2 = 0), we get 

(2.13) C{dg)Uf'' = -di (903d3 djip 3 {dAj{ 0 ) ■ n) r 3 - 03 dg^oi djcpi (dAj(O) • r 3 ) ri. 

In particular, the spectrum of is a subset of the integers a G that satisfy 0:2 = 0 and ai 0(3 7 ^ 0, so 
has zero mean when integrated with respect to ( 0 i, 03 ). 

Equation (2.6) (b) also shows that the component Ui that carries the 02-oscillations of Ui satisfies 
L{d(p 2 ) dg^Ul = 0, so that Uf can be written as = T 2 {t,x, 02 )r 2 for an unknown scalar function T 2 of 
zero mean with respect to 62 - 

Let us now consider Equation (2.6) (c). Since Uq only has oscillations in 0i and O 3 , and since there is 
no resonance among the phases, none of the terms Ni{Uq,Uq), M 2 {Uo,Uo,Uo) has oscillations in 02 only. 
Looking also closely at each term in M{Uq,Ui),M{Ui,Uq), we find that both expressions have zero mean 
with respect to ( 0 i, 03 ), because the only way to generate a 02 -oscillation would be to have a nonzero 
mode of the form (ai,a 2 , 0 ) or { 0 , 02 , 03 ) with 02 / 0 in Uf^, but there is no such mode according to 

(2.13) . We thus derive the outgoing transport equation 

dm + V2 • Vm = 0 , 


from which we get T 2 = 0 . 
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Step 4: computation of the nonpolarized components of U.\,1A\^ and compatibility condition. 
At this stage, we know that the first corrector Ux reads 


Ux =Ux{t,x,dx-,h) = U.x{t^x) + Ul{t,x,ex) + Ul{t,x,9z) + x,^s), 

with tlf^ determined by (2.13). Moreover, the nonpolarized part of is obtained by considering 
Equation (2.6) (b) and retaining only the @ 1^3 Fourier modes. We get 

E(d(/9m) dg^lAx — ^{9') {Ofxi fm) ^ra dji^rn (dAj (0) • , Ul — 1,3, 

SO (/ — m = 1,3, is the only zero mean function that satisfies 

(2.14) (/ - Pm) deJAY" = -Rm L{d) (am Xm) - CTm de^am djipm Rm (dAj(O) ■Xm)Xm, m = 1, 3 . 

We now consider the boundary condition (2.8) (b), which we rewrite equivalently as: 

BU.x\xd=o + B Pi U^=o,ei=6»o + B P‘iUl\xj^=ofi^=eQ 
+B (I - Px)lAx\xd=o,01=00 + B (I - P3)Uf\x^=o,d3=eo + BUf‘^\x^=ofii=03=0o 

(d26(0) • (e, e)) = G(t, y, 9o). 

We differentiate the latter equation with respect to 9o and apply the row vector b, so that the first line 
vanishes. We are left with 

bB (I - Px) (deM)\^d=o, 9 i= 0 o + bB (I - P 3 ) {d 0 ,Uf)\x,=o, 93 = 0 o+kBd 0 ,(Uf%^=o, 0 i= 03 = 0 o) 

+ ^k{<i^b(0) ■ (e,e)) 90 o(a^) = bd0^G. 

The first two terms in the first row are computed by using (2.14), and [CGIO, Proposition 3.5]. We get 

(2.15) vde^(a^) — Xi^o-pO- + kBd 0 Q{Ux'^\xd=o, 01 = 03 = 00 ) = bde^G, 
where the constant v and the vector field Xlop are defined by 

(2.16) 

:= ^b(d%(0) ■ (e,e)) - ^bBRxdjipx (dAj(O) - 61)61 - ^bBR^dj^p^ (dAj(O) - 63 ) 63 , 

(2.17) 

d-l 

Xlop := bB (Rx 61 + R3 63) dt + ^bB (Rx Aj{0) 61 + R3 Aj{0) 63) dj = i ( 5 ^<T(r,^) dt + dr^^a{T,y) dj), 

i=i 

with L a nonzero real constant, and a defined in Assumption 1.6. It is also shown in [CGIO, Proposition 3.5] 
that the partial derivative dra(T,ri) does not vanish, so that, up to a nonzero constant, Xuop = dt + 'w■'Vy 
for some vector w G (which represents the group velocity associated with the characteristic set of 

the Lopatinskii determinant). 

Step 5: computation of Uf^ and conclusion. 
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The final step in the analysis is to compute the derivative 61 = 63 = 60 ) arising in (2.15) in 

terms of the amplitude a. Restricting (2.13) to the boundary {xd = 0} gives 

C{de)U^%^=Q = - a(t,y,6'i) {dgga){t,y,e 3 ) djtps (dRj(O) • 61)63 
- a{t,y,e 3 ) {dgga){t,y,ei)dj<fi (dAj(O) • 63)61 . 

Let us expand a in Fourier series with respect to 6 q (recall that a has mean zero): 

a{t,y,eo)= 

keZ* 


Then the Fourier series of reads 


with 


(2.18) 


L/r(Lx,0i,03)= Ufc,,fc 3 (La;)e 2 *-("i®i+" 3 e 3 )/ 0 ^ 

fcijfcsGZ* 


Uki,k3{t, y, 0) = -L{ki dipi + k3 dvJs) ^ (fci + ks Es^i ), 

El,3 := djipi (d^j(O) • 63) 61 , £’3,1 := dj(p3 (d^j(O) • 61) 63 . 


Plugging the latter expression in (2.15), we end up with the evolution equation that governs the leading 
amplitude o on the boundary: 


(2.19) 


) ^Lop<^ “t“ 9^0 Qper [^5 


with 


Qper [o), tt] i — 


E 

fcGZ 




6£L(A:id(/?i + A:3d(/?3) ^ (/ci £1,3 + A:3 £3,1) 0^,3 


ki+k 3 =k, 

ki k 3 j !^0 


^2 i TT A: 60 /0 


/ 


Equation (2.19) is a closed equation for the leading amplitude a on the boundary. It involves the vector 
field Xlop associated with a characteristic of the Lopatinskii determinant, a Burgers term dg^a^ and a new 
quadratic nonlinearity dg^Qper{<i-, a]- The operator takes the form of a bilinear Fourier multiplier. Its 
above expression may be simplified a little bit by computing the decomposition of L(A:i d(^i + k 3 dip 3 )~^ 
on the basis ri,r 2 ,r 3 , and by recalling the property hBri = bBr 3 = 0 (so only the component of 
L{ki d(pi + k 3 dnp 3 )~^ on r 2 matters). We obtain: 


( 2 . 20 ) 


Qper[a,a] ■= -bBr2 


E 

k£Z 


ki h £1,3 + ^3 h £3, 


1 


y _ 

k ^^ k ^= k , ^-1 “ -2) + ^3 - ^2) 


^ki 


\ k^k^jtO 


,^2 i TT fc 60/0 


/ 


Anticipating slightly our discussion in Section 3, well-posedness of (2.19) will be linked to arithmetic 
properties of the phases t-pm^ and this is one reason for the small divisor condition in Assumption 1.9. This 
is in sharp contrast with the theory of weakly nonlinear geometric optics for both the Cauchy problem 
(see [HMR 86 , JMR93, JMR95] and references therein) and for uniformly stable boundary value problems 
(see [Wil99, Wil02, CGWll]), where arithmetic properties of the phases do not enter the discussion on 
the leading profile for the high frequency limit. 
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2.3 Extension to strictly hyperbolic systems of size N 

The above derivation of Equation (2.19) can be extended without any difficulty to the case of a hyperbolic 
system of size N provided that Assumption 1.9 is satisfied. In that case, the number of phases equals N. 

Steps 1 and 2 of the above analysis extend almost word for word, to the price of changing some 
notation. Namely, the first relations of the WKB cascade (2.6), (2.8) shows that the leading amplitude 
Uq reads 

N 

... ,6 n) =lio{t,x) + ^ am{t,x, 0 m)rm , 

m=l 

with unknown scalar functions depending on a single periodic variable 9^, and of mean zero with respect 
to this variable. The quadratic expression A4{Uq,Uo) still has zero mean with respect to (0i,..., 0Ar) so 
the nonoscillating part Wq satisfies (2.10) and (2.11), and therefore vanishes. Furthermore, each function 
am satisfies the Burgers equation (2.12), which reduces the leading amplitude Uq to 

(2.21) Lioit, X, 01, ... , 07v) — ^ ^ *^m(b 9m) j 

m£X 

where I denotes the set of incoming phases. The boundary condition (2.8) (a) then yields 

\fmeZ, amit,y, 0, 9o) Vm = a{t, y, 9o) , 

for a single unknown scalar function a of zero mean with respect to its last argument 9o (e = Xlm^x 
spans the vector space E®(r, rj) D Ker B). 

Step 3 of the above discussion is unchanged, showing that in the first corrector Ui, each profile 
vanishes when the index m corresponds to an outgoing phase. The noncharacteristic part is obtained 
by using the relation 

£idg)Ur = - E ^&7n2^'^2 ^j^m2 (®) ‘ ^mi) ^m2 ^m2 ^mi djfm\ (dAj(O) • Tm 2 ) • 

m\<m 2 

mi ,m2 GX 

which is the analogue of (2.13). 

Step 4 is also unchanged because lAi has no outgoing mode, and when m corresponds to an incom¬ 
ing phase, (/ — Pm)d 0 ^U^ is given by (2.14). Eventually, the boundary condition (2.8)(b) gives the 
compatibility condition 

( 2 . 22 ) V dg^{a^) - Xi^opa + dgoQper[a,a] = bdg^G , 

with 

(2.23) V := ^b{d^b{0) ■ {e,e)) - ^ ^ bB Rmdjipm idAj{0) ■ Cm) em , 

7n£X 

d-l 

^Lop • — ^ ^ ^ B Rm EE bBRm Aj (0) Cm dj = i {draiz, y) dt + dr^^a{T, y) dj ), 

m£X j=l m^X 

where t is a nonzero real constant and the function a is defined in Assumption 1.6. (Again, [CGIO, 
Proposition 3.5] shows that the partial derivative dra{T, y) does not vanish.) The new expression of the 
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bilinear Fourier multiplier Qper reads: 


(2.24) Qper[o,a] := -'^hBr^ ^ 


mec> mi<m2 

mi ,m2 GX 


/ 




E E 



^km2 ® 


^2i7r kOo/S 


ks^ kmi~\~km2 — k, 

\ ^mi 


kmi (^mi ^m) + ^m2 


/ 


with 


(2.25) 


V mi, m 2 € X, FI, 


•— dj^mi (d-^j(O) ■ ('m2) (mi • 


The definition (2.24) reduces to (2.20) when N = 3 and Assumption 2.3 is satisfied. 

2.4 The case with a single incoming phase 

In this short paragraph, we explain why the computations in [AM87, WY14] lead to the standard Burgers 
equation for determining the leading amplitude, and does not incorporate any quadratic nonlinearity of 
the form (2.20) we have found under Assumption 2.3. 

The vortex sheets problem considered in [AM87] and the analogous one in [WY14] differ from the 
framework that we consider here by the fact that the (free) boundary is characteristic. Nevertheless, 
one can reproduce a similar normal modes analysis for trying to detect violent or neutral instabilities. 
The two-dimensional supersonic regime considered in [AM87] precludes violent instabilities, but a similar 
situation to the one encoded in Assumption 1.6 occurs®. The situation in [WY14] for three dimensional 
steady flows is similar, and the corresponding Lopatinskii determinant is computed in [WY13]. 

The two-dimensional Euler equations form a system of three equations {N = 3), but due to the 
characteristic boundary (the corresponding matrix A,i(0) has a kernel of dimension 1), the number of 
phases ipm on either side of the vortex sheet equals 2. One of them is incoming, and the other is outgoing. 
In such a situation, there are too few incoming phases to create a nontrivial component Uf^ for the first 
corrector Ui, so that the bilinear Fourier multiplier Qper vanishes. Though our argument is somehow 
formal, the reader can follow the computations in [AM87] or in [WY14] and check that they follow the 
exact same procedure that we have described in our general framework. 


3 Analysis of the leading amplitude equation 


Our goal in this section is to prove a well-posedness result for the leading amplitude equation (2.19). Up 
to dividing by nonzero constants, and using the shorter notation 6 instead of 00 , the equation takes the 
form 


(3.1) 


dta-\-w ■ Vytt + cadga pdgQpeila, a] = g , 


®The reader will find in [CS04] a detailed analysis of the roots of the associated Lopatinskii determinant, showing that 
they are located in the hyperbolic region and simple. 
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where w is a fixed vector in are real constants, and Qper is a bilinear Fourier multiplier with 

respect to the periodic variable 0 \ 


(3.2) 


Qper o] • — 


E 


E 


ki+k 3 =k, 


h h -£^1,3 + fcs h ^3,1 

fei (wi - W2) + kz (^3 - W2) 


^ki ^^3 


^2 ZTT /c 0/© 




The source term g in (3.1) belongs to H^°°{{—oo,TQ\t x x (]R/(0Z))^), Tq > 0, and vanishes for 

t < 0. Furthermore, it has mean zero with respect to the variable 6 . Recall that in (3.2), denotes the 
k-lh Fourier coefficient of a with respect to 9 (which is a function of {t,y)). 

Recall that for strictly hyperbolic systems of size N, (2.24) should be substituted for (3.2) in the 
definition of Qpei, while in the particular case p = 1 (one single incoming phase), (3.1) reduces to the 
standard Burgers equation for which our main well-posedness result, Theorem 3.4 below, is well-known. 
For simplicity, we thus encompass all cases by studying (3.1), (3.2) and leave to the reader the very minor 
modifications required for the general case. 


3.1 Preliminary reductions 

We first introduce the nonzero parameters: 

53 := 

that satisfy ^3 = 1 + (5i, and we observe that Qper[a, a] in (3.2) can be written as 

r„ „i _ 0^2 -E'i,3 Ttp ta „ & h Es,! 


2 vr (^3 — W;t) 
where the bilinear operator Fper is defined by: 


/ 


(3.3) 


IP'per(^)^) •— 'y ^ 
keZ 


2 TT (^3 — ) 


E 


ki+k 3 =k, 


ki + ks S 3 


^2ITT A: 6/0 




The bilinear operator Fper satisfies the following two properties: 


(3.4) 


(Differentiation) d 0 (Fper(u, v)) = Fper(d 0 U, v) + Fper(u, Ogv) 


(3.5) (Integration by parts) Fper(tt, dev) = --^^u u - ^ Fper(<96)M, u), if uq = uq = 0 . 

0 03 03 


Using the properties (3.4), (3.5), we can rewrite equation (3.1) as 


(3.6) dta + w • VyO + cadga + /rFper(den) dgo) = g , 

with new (harmless) constants c, g for which we keep the same notation. 0 ur goal is to solve equation 

(3.6) , that is equivalent to (3.1), by a standard fixed point argument. The main ingredient in the proof 
is to show that the nonlinear term Fper(d 0 fl) dga) acts as a semilinear term in the scale of Sobolev spaces. 


20 













3.2 Tame boundedness of the bilinear operator Fper 

The operator Fper is not symmetric but changing the roles of di and 63 , the roles of the first and second 
argument of Fper in the estimates below can be exchanged. This will be used at one point in the analysis 
below. We let := x (1R/0Z)) denote the standard Sobolev space of index € N. The norm 

is denoted || • \\h''- Functions are assumed to take real values. In the proof of Theorem 3.1 below, we 
shall also make use of fractional Sobolev spaces on the torus M/0Z or on the whole space These 

are dehned by means of the Fourier transform, see, e.g., [BGS07, BCDll]. Our main boundedness result 
for the operator Fper reads as follows. 

Theorem 3.1. There exists an integer vq > 1 + d/2 such that, for all v > vq, there exists a constant Cy 
satisfying 

(3.7) Vu,v G ff", ||Fper(d 0 U, < Cy (||u||h-o ||u||h- + ||w||h- Ibib-o) . 

Estimate (3.7) is tame because the integer ng is fixed and the right hand side of the inequality depends 
linearly on the norms ||u||//i', This will be used in the proof of Theorem 3.4 below for propagating 

the regularity of the initial condition for (3.6) on a fixed time interval. 

Proof. We first observe that, provided that Fper(u, u) makes sense, then Fper(u, u) takes real values. This 
simply follows from observing that 


(Fper('W, u))_^ — (Fper('U., u))/j , 

provided that u and v take real values. We are now going to prove a convenient new formulation of 
Assumption 1.9. 

Lemma 3.2. There exist a constant C > 0 and a real number 70 > 0 such that 

Vfci,/c 3 gZ\{ 0 }, ^^3 J 3 I < C mm{\kip° ,\k3p°). 

Proof of Lemma 3.2. Since in our framework, we have X = {1, 3} and O = {2}, we can apply Assumption 
1.9 to any (fci, 0, k^) G I? with ki fcs 7 ^ 0. We compute 


L(d(A:i ipi + k3 (ps)) r2 = (ki (w^ - ^ 2 ) + k^ (wg - ^ 2 )) Ad(0) r2 , 


and the quantity ki (w^ — ^ 2 ) + A :3 (wg — W 2 ) cannot vanish for otherwise there would be a nonzero vector 
in the kernel of L(d(ki ipi + /cg yjg)). We thus derive the bound 


\ki (wi 


1 

(^2)+ ks (Wg 


kl2)l 


< C\\L{d{ki(pi + k3(p3)) ^11 , 


for a suitable constant C that does not depend on /ci, fcg. The norm oi L(d{ki ipi + k^ (p 3 ))~^ is estimated 
by combining the lower bound given in Assumption 1.9 for the determinant, and an obvious polynomial 
bound for the transpose of the comatrix. We have thus shown that there exists a constant C > 0 and a 
real parameter 70 (which can be chosen nonnegative without loss of generality), that do not depend on 
/ci, fcg, such that 

nn-7X177- g<C'|(fci,A:3)r . 

1^1 (kli ~ kl2) + ^3 (kls “ 1 X 2)1 
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Up to changing the constant C, we can rephrase this estimate in terms of the rescaled parameters (5i, 


(3.8) 


1 


\ki 


<C\{kuhW% 


and it only remains to substitute the minimum of |/ci|, j/caj for the norm |(A:i, ^ 3 )! in the right hand side 
of (3.8). 

There are two cases. Either \ki6i + k^d^l > |5i| > 0, and in that case, it is sufficient to choose 
C > l/|5i| (we use 70 > 0 ). Or + A; 3 ( 53 | < |(5i|, and we have 

|A:3| < ^ \ki <5i + fca <^3| + ^ 1^1 5i| < 2 ^ Ifcil, 
because ki is nonzero. Up to choosing a new constant C, (3.8) reduces to 

^ <C\ki\^\ 


\ki + fcaJal 

and we can prove the analogous estimate with fca instead of ki by the same arguments. This completes 
the proof of Lemma 3.2 □ 

The proof of Theorem 3.1 relies on the following straightforward extension of [RR82, Lemma 1.2.2]. 
The proof of Lemma 3.3 is exactly the same as that of [RR82, Lemma 1.2.2], and is therefore omitted. 

X Z X X Z —)■ C 6e a locally integrable measurable function such that, 


Lemma 3.3. Let 

either 


or 


Then the map 


dr] <+00, 


<+00- 


{f,g )'—^ / '^^{^,k,r],£)f{^-ri,k-£)g{7],i)dr], 

.htd-l 


i&'L 


is bounded on L^(M'^ ^ x Z) x L^(M'^ ^ x Z) with values in L^(M'^ ^ x Z). 

To prove boundedness of the bilinear operator Fper(i96C) de-), we shall apply Lemma 3.3 in the Fourier 
variables. More precisely, for functions u,v in the Schwartz space x (M/0Z)), there holds'^: 


Cfc(Fper(<9eu,<9eu))(0 = C 


St 




{k-£)£ 


ik-£) 61+£63 


Ck-i{u){f - rf) ce{v){r]) dry. 


e£Z,e^{o,k} 

Omitting from now on the constant multiplicative factor, we consider the symbol 


K{k,£) := 


{k-£)£/{ik-£) 6 i+£ 53 ) d£^{d,k}, 

0 otherwise. 


'Here we use the notation Ck for the fc-th Fourier coefficient with respect to the variable 6, and the ”hat” notation for the 
Fourier transform with respect to the variable y. 


22 






We wish to bound the H'^ norm; 

[ E dev))m^ , 

for 1 / G N large enough ((•) stands as usual for the Japanese bracket). 

Given the parameter 70 > 0 in Lemma 3.2, we hx an integer t'g > 7o + 2 + d/2. We consider two 
functions Xi, X 2 on X such that xi + X 2 = Ij and 

Xi{^,k,r],i) = 0 if {{T],i)) > (2/3) {{C,k)), 

X 2 (C, k, 7], £) = 0 if ((77, £)) < (1/3) ((^, k )). 


We first consider the quantity 

(3.9) [ xi{C,k,r],e) {{^,k)y K{k,e)ck-eiu){C - r])ceiv){r])di], 

which we rewrite as 

[ ^ xi{C,k,v,£){{C,k)rK{k,e) 

' {i^-v,k-£)r{{v,£)yo 


-r],k- £)Y Ck-iiu) (? - 7 ?)) (((7, ^)y° ci{v) ( 7 )) dT?. 


On the support of xi, there holds 


{{^-V,k-i))>{{^,k))-{{riY))>^{{^,k)). 


and therefore 


E 

iez 


Xi{^,k,vJ) {{Yk)Y K{k,£) 


{{^-r],k-£)Y {{7 ],£)Yo 


*' s c I y: 


\K{k,£Y 


We now use Lemma 3.2 to derive the bound 


from which we get 


{k-£)£ 

1 


{k-£) J1+M3 

l< 5 i| 

{k-£) J1+M3 


sup 


[ E 

■ixZ JK'* 1 


Xi{^,k,vY){i^,k)mk,i) 


- Tl,k - £)Y {{7 ],£)Yo 


dr] <C 


< C |fP“+^, 

|£|2(70+2) 


d77. 




iivY)) 


2 UO 


dr] < + 00 , 


d^ E G 11 77 11 11 "^11/7^0 


thanks to our choice of vq. Applying Lemma 3.3 to the quantity in (3.9), we obtain 

2 

/ E / '^Xii^,k,riY){iC,k)Y't^{k,£)ck-e{u)iC-v)ce{v){r])dri 

Similar arguments yield the bound 

/ E / '^X2{C,k,vY){{Yk)YK{k,£)ck-i{u){C-r])ci{v){r])dT] 

and the combination of the two previous inequalities gives the expected estimate 

/ E ((?>^))^^ |Cfe(Fper((9(7«, l96)7;))(^)|^d^ < C (IIuIIh'^o + WuWh’' lbll/7‘"o) 


d^ < C llttll^i^o ll^^l 


H'' ) 


□ 
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3.3 The iteration scheme 


In view of the boundedness property proved in Theorem 3.1, Equation (3.6) is a semilinear perturbation 
of the Burgers equation (the transport term w • Vy is harmless), and it is absolutely not surprising that 
we can solve (3.6) by using the standard energy method with a fixed point iteration. This well-posedness 
result can be summarized in the following Theorem. 

Theorem 3.4. Let uq be defined as in Theorem 3.1, and let v > vq. Then for all R > 0, there exists a 
time T > 0 such that for all data oq G x (M/0Z)) satisfying ||ao||_ff‘"o < R, there exists a unique 

solution a G C{[0,T]; H'^) to the Cauchy problem: 

{ dto + w • VyO + cadga + fj,¥per{dga, dga) = 0 , 

a\t=o = do ■ 

In particular, if oq G x (M/0Z)), then a G C{[0,T]; where the time T > 0 only depends 

on llaoll/f'^o • 

Proof. The proof follows the standard strategy for quasilinear symmetric systems, see for instance [BGS07, 
chapter 10] or [Tay97, chapter 16], and we solve the Cauchy problem by the iteration scheme 

J+ w • + ca” + fj.¥pe,.{dga‘^, dgo^) = 0 , 

= ao,n+i, 

where (oo,n) is a sequence of, say, Schwartz functions that converges towards oq in H'', and the scheme 
is initialized with the choice = ao,o- Given the radius R for the ball in we can choose some time 
T > 0, that only depends on R and z/, such that the sequence (a"') is bounded in C([0, T]; 77^). The 
uniform bound in C([0,T];77^) is proved by following the exact same ingredients as in the case of the 
Burgers equation. Contraction in C([0,r];L^) is obtained by computing the equation for the difference 
j.n +1 ._ g^n+i _ which reads 

(3.10) + w • + ca" dgr^+^ = -cr^ dgol^ - pL¥p^fidgr^, dgoJ^) - pL¥p^fidgaI^-^,dgr^). 


The error terms on the right hand-side are written as 

Fper(5er", dgal^) = -|^ r" dga^ - | Fper(r", a|ea"), 

¥p,fidga^-\der^) = -^r'^dga^-^ _ ^ Fper(5|,a"-\r"), 

fe) 03 O 3 

where we have used (3.5). 

The final ingredient in the proof is a continuity estimate of the form 


(3.11) |lFper(u,u)lli 2 < C min [\\u\\h<' 0-2 I^IIlz IblliT-o-a) , 

which we now prove for completeness. We apply the Fubini and Parseval-Bessel Theorems to obtain 

2 


llFper(?^, v)\\i ^2 



E 

ki+k3=k, 


1 

ki (5i -7 ks 63 


'^ki 


dy, 
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and then apply the -k P ^ 1“^ continuity estimate to derive 

\Wperiu,v)\\l 2 <C j l«fc| I 

^ VfcGZ / fcGZ 

We then apply the Cauchy-Schwarz inequality and derive the estimate 

||lF'per('U, u) 11^2 < C* f sup ||'u(y,') ll/f 7 o+l(R/© Z) ) 11 '^IIl 2 ) 

which yields (3.11) because the integer vq in Theorem 3.1 can be chosen larger than [d — l)/2 + 70 + 3. 
(The ’’symmetric” estimate is obtained by exchanging the roles of u and v.) 

At this stage, we multiply Equation (3.10) by and perform integration by parts to derive 

sup ||r’^+ii 2 < ||r’"+^| 4 =o|li 2 +C'or sup \\r^+^ \\l2 + Cq T sup \K+^\\l 2 sup ||r"||L2 , 
te[o,T] te[o,r] te[o,T] te[o,T] 

where the constant Cq is independent of n and follows from the uniform bound for sup^gjo.T] 

By classical interpolation arguments, (a”) converges towards a weakly in L°°{[0,T]‘, H'^) and strongly in 
< u. Continuity of a with values in is recovered by the standard arguments, see, e.g., 
[Tay97, Proposition 1.4]. 

If 1 / > z^o, it remains to show that the time T only depends on the norm ||ao||//‘'o, and this is where the 
tame estimate of Theorem 3.1 enters the game. More precisely, we follow the same strategy as in [Tay97, 
Corollary 1.6], and show that the Tf^^-norm of the solution a satisfies a differential inequality of the form 

4hA!E<c„(iio(t)ii7.)i|n(t)iii,., 

where Cy is an increasing function of its argument. In particular, boundedness of a{t) in 77^° on an 
interval [0,T'), T' > 0, implies a unique extension of the solution a G C([0, T'); 77^) beyond the time T\ 
which means that the time T of existence for a only depends on UaoljH^'o. □ 

3.4 Construction of the leading profile 

Theorem 3.4 is the cornerstone of the construction of the leading profile T/q. Solvability of (2.22) for a is 
summarized in the following result. Recall that the smoothness assumption for G was made in Theorem 
1 . 10 . 

Corollary 3.5. There exists T > 0, and a ^ C°°{{—oo,T]] ^x(M/0Z))) solution to (2.22) with 

fl|t<o = 0- Furthermore, a has mean value zero with respect to the variable 9. 

Proof. Equation (2.22) is easier to solve than the pure Cauchy problem in Theorem 3.4 because we can 
apply Duhamel’s formula starting from the initial condition oq = 0. From the assumption of Theorem 
1.10, we have G G C°°((— 00 ,Tq]; 77+°°(M'^“^ x (M/0Z))) with Tq > 0 and Gjt<o = 0, so we can find 
0 < T < To and a G C((— 00 , T]; 77+°°(M'’*“^ x (M/0Z))) solution to (2.22) with ajf<o = 0. Here the time 
T depends on a fixed norm of G. Then Equation (2.22) yields a G C°°((— 00 , T]; 77+°°(R'^“^ x (R/0Z))) 
by the standard bootstrap argument and smoothness of G. 

For every fixed y, the mean value 

1 7® 

a(t,y)-=-Q a{t,y,6)de, 
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satisfies the homogeneous transport equation 


dta + w • VyO = 0 


with zero initial condition, and therefore vanishes. 


□ 


After constructing a on the boundary, we can achieve the construction of the leading profile Uq in the 


whole domain. 

Corollary 3.6. Up to retricting T > 0 in Corollary 3.5, for all m & I, there exists a unique solution 
am G C°°{{—oo,T]; x (M/0Z))) to (2.12) with am\t<o = 0 and am\xd=o = where the real 

number tm is defined by Cm = tm^m- Furthermore, each am has mean value zero with respect to the 
variable 9m- 

The result follows from solving the boundary value problem for the Burgers equation (2.12) with 
prescribed Dirichlet boundary condition on {xd = 0}. This is a (very!) particular case of a quasilinear 
hyperbolic system with strictly dissipative boundary conditions for which well-posedness follows from the 
classical theory, see, e.g., [BGS07]. 


The leading profile Uq is then given by (2.21) and belongs to C°°((—oo, T]; x (]R/0Z)-^)). 

Furthermore, its spectrum with respect to the periodic variables (0i,..., Op^) is included in the set 


:= {a € jMm £ O, am = O} 


(3.12) 


as claimed in Theorem 1.10. 


4 Proof of Theorem 1.10 


4.1 The WKB cascade 


In this paragraph, we give a more detailed version of (2.6)-(2.8). We plug again the ansatz (2.5) in (1.1) 
and derive the set of equations (4.1)-(4.3) below. We recall that the operators C,A4 are defined in (2.7), 
while L{d) is defined in (1.4). Then the WKB cascade in the interior reads: 


(a) C{de)Uo = 0, 


(4.1) (b) £idg)lfi + L{d)Uo + MiUo,Uo) = 0, 

(c) £{d0)Uk+2 F L{d)Uk+i + FA{Uo,Uk+i) + M{ldk+i,hio) +Fk = 0, k>0 


with 


(4.2) 






• • 5 ) 1 


£=1 K,l-\ - \-K£ = U—£ 
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Observe that (4.1)(c) coincides with (2.6)(c) for k = 0. Furthermore, each matrix AJ, i' > 1, only depends 
onUo,... and therefore each source term F^, k > 0, only depends on Uq, ... 

The set of boundary conditions for (4.1) reads (recall B = d5(0)): 



(a) 

B IAq — 0 , 


(4.3) 

(b) 

BUi + ]^d?biQ,)- 

{Uo,Uo) = Git,y,9o), 


(c) 

BUk+2 + d%it)) ■ 

{lAQ,lAk+i) + ^k = 0 , A: > 0 , 

with 




(4.4) 

V/c > 0, 

k+1, 

E d'6(0)-(Zi.i,...,77,, 

K,i-\ - = 


In (4.3) and (4.4), all functions on the left hand side are evaluated at = O,0i = ■ ■ ■ = On = Oq- The 
source term in (4.3) (c) only depends on Uq, ... ,Uk- 

We are looking for a sequence of profiles {Uk)ken that satisfies (4.1)-(4.3), and Uk\t<o = 0 for all k. 

4.2 Construction of correctors 

Some notation will be useful in the arguments below. For any function / that depends on , 0^), 

with 0-periodicity with respect to each 6^, we decompose / as 

N 

f = lit, x)+Y^ f^it, X, 0^) + /“(t, x, 01, ..., 0iv) , 

m=l 

where / stands for the mean value of / on the torus (M/0Z)'^, each incorporates the 0m-niodes 
of / (in particular, the spectrum of is included in and fm has mean zero with respect to 

9m), and the spectrum of is included in \ Here, the spectrum only refers to the Fourier 

decomposition of / with respect to (0i,... , 0Ar). The mappings / i—>■ and / i—are continuous 

on C°°((—oo, T]; x (M/0Z)'^)). Furthermore, if the spectrum of / is included in Z^, then 

belongs to the space of profiles defined in Lemma 4.1 below. 

The following observation is well-known in the theory of geometric optics, see for instance [JMR93, 
Wil99], and relies on Assumption 1.9. 

Lemma 4.1. The operator Cidg) is a bounded isomorphism from into itself, where 

P“ := |/ e C°°((-oo,r];F+°°(M^ x (M/0 Z)^)) / Spectrum (/) C Z^ . 

Indeed, for a € Z^ \ Z'^d^ matrix L(d(a • <h)) is invertible and the norm of its inverse is bounded 

polynomially in \a\ (the degree of the polynomial being fixed). We shall feel free to write C{de)~^ /“ 
when is an element of P'^^. 

Unsurprisingly, the construction of the sequence iUk)k&i is based on an induction process. We formu¬ 
late our induction assumption. 

(H„) There exist profiles Uq, ... ,Un in C°°((—oo,T]; x (M/0Z)^)) that vanish for t < 0, whose 
( 01 ,..., 0iv)-spectrum is included in Z^, and that satisfies 
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• (4.1) (a) and (4.3) (a) if n = 0, 

• (4.1)(a)-(b) and (4.3)(a)-(b) if n = 1, 

• (4.1)(a)-(b), (4.3)(a)-(b), (4.1)(c) and (4.3)(c) up to order n — 2 if n > 2, 

• Compatibility condition in the interior: 

(4.5) VmGX, = and ^ = 0, 

• Compatibility condition on the boundary: 

(4.6) 6(- 

vri£X 

+ dg^Gnit, y, 0 q)J = 0 , 


where in (4.5) and (4.6), we have set: 


(4.7) 

Fn- = 

(4.8) 

Gn : = 


L{d)Uo + M{Uo,Uo), ifn = 0, 

Lid)Un + MiUo,lln)+Mil{nM)+fn-l, if n > 1 , 

i d'^b{0) ■ {Uo,Uo)\x^=o,di=-=gj^= 0 o - G{t, y, 6*0), if n = 0 , 
d'^b{0) ■ {Uo,Un)\xa=o,ei=-=eN=eo + ^^n-i, if n > 1. 


Recall that and Gk are dehned in (4.2) and (4.4), so the above source terms Fn,Gn only depend on 
Uo,...,Un- Several properties of these source terms are made precise below which, in particular, will 
justify why we can apply the operator C{dg)~^ to 

The reader can verify that our construction of the leading profile Uq in Sections 2 and 3 proves that 
(Ho) holds for some T > 0. In that case, (4.6) reduces to (2.22), which was solved in Section 3. The 
compatibility conditions (4.5) in the interior give the decoupled Burgers equations (2.12) for each incoming 
amplitude am- 

Our goal is to show that (H„) implies (H„,+i) with the same time T > 0, which will imply that there 
exists a sequence of prohles {Uk)keN in C°°{{—c>o,T]-, x (M/OZ)-^)) that satisfies (4.1)-(4.3), and 

b{k\t<o = 0 for all k. We decompose the analysis in several steps, as in Section 2, assuming from now on 
that (H„) holds for some integer n € N. 

Step 1: properties of Gn, and definition of (/ — Pm)^^i- 

From assumption (H„), the profiles Uq, ... ,Un belong to C°°((—00, T]; x (M/0 Z)^)), vanish 

for t < 0, and their spectrum is a subset of Z^. The space of such functions is an algebra, and therefore, 
we can verify from (4.2) and (4.7) that Fn belongs to C°°{{—oo,T]; x (M/GZ)^)), vanishes for 

t < 0, and its spectrum is a subset of Zj . Consequently, F^'^ belongs to the space of profiles defined 
in Lemma 4.1. 

Our goal is to construct a prohle Un+i that satishes 
(4.9) C{dg)Un+l + Fn = 0 . 

We first observe that if a € Z-^ is a noncharacteristic mode, that is, a € Z-^ \ Z^’^, and if moreover a 
has one nonzero coordinate am with m (z O, then the a-Fourier coefficient of Un+i vanishes. This implies 
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that any solution to (4.9) has its spectrum included in Zj . We thus define := —C{do)~^ so 

Un+i belongs to P“ and vanishes for t < 0. 

For all m = 1,..., A^, we define (I — Pm) as the unique mean zero solution to 

In particular, we can write '^rn for m € O since has no outgoing mode. Due to the 

compatibility condition (4.5) for m G Z, we have 


Vm = 1,..., iV, Lidipm) (/ - Pm) deX+1 +K = 0, 


which means that the components of Un+i that we have already defined satisfy 


C{de) 


N 


\m=l 


{I-Pm)K+l+U, 


'‘n+1 


+ Zn — 0 , 


because Fn has zero mean value. It is clear that the nonpolarized components (/ — belong to 

C°°((— CX 3 , T]; X (M/GZ)-^)) and vanish for t < 0, because the F^’s do so. 


Step 2: Un+i has no outgoing mode. 

Let m £ O. We focus on (4.1) (c) for k = n. Since has no outgoing mode, the profile Un+i must 
necessarily satisfy 

immi ^n+1^m) P Im {M{U^Mn^l) + M{Un^x,U^))'^ = 0 • 

However, the leading profile is given in (2.21), and since the only noncharacteristic modes of belong 
to Zj , we observe that it is not possible to generate a 6*m-iiiode in M.{UQ^Un+i) nor in M.{U^_^i^Uq). This 
means that the amplitude cr((Yi satisfies the outgoing transport equation 

dtfJn+l + Vm • ^x(Tn+l = ^ > 

and therefore vanishes. 

Since the profile Un+i has no outgoing mode, it reads 

Un+l = ^ ^ Z^r+l + E ^n+l + Kh , 

m^X mGX 

and it only remains to determine the mean value Un+i and the polarized components Pml^^+i, m £ I. 
Let us observe that such components have no influence on the fulfillment of (4.9), no matter how we 
define them because they belong to the kernel of €{ 80 ). Hence we shall no longer focus on (4.9), but 
rather on the boundary conditions for Un+i and the interior compatibility condition at the next order. 


Step 3; determining 

Let us first derive the interior equation for Un+i- We introduce the notation 
m £X. We consider (4.1) (c) for k = n, and take its mean value on the torus, observing first that both 
terms M.{Uq,U^‘^i) and M.{U^‘^i,Uq) have zero mean value. Hence we derive the equation 

rS 

L{d) Un+l + y] / O'm -PT ^iZm (d^i(O) • Tm) Tm 

+ y / CF'^+l^0^am-^^j^Pm{<^Aj{{)) ■rm)rm +Xn{t,x) 
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where the source term is defined by: 


7-^ := F^+ ^ (/ - • 

m^X Tn£X 

We observe that each sum of integrals 

vanishes, and therefore Un+i must satisfy the system 

L{d)Un+l=-Tn, 


in (—oo,T] X The source term Tn belongs to C°°((—oo, T];and vanishes for t < 0. 

The boundary conditions for Un+i are obtained by taking the mean value of (4.3) (b) if n = 0 or 
(4.3) (c) for /c = n — 1 if n > 1. In any case we find that Un+i\xd=o must satisfy 

B lAn+l |xd=0 T B \xd=0,6i=---=9M=9o ) 4“ Gn jt, t/) = 0 . 

Since already been determined, we can apply the well-posedness result of [Cou05], supplemented 

with the regularity result in [MSll], and construct a solution Un+i G C°°((—oo, T]; ff+°°(M!^)) to the above 
eqnations (in the interior and on the boundary). 

Step 4; determining Part I. 

We keep the notation for m G Z. The evolution of is obtained by imposing 

the compatibility condition: 

£™L(a)i^;r+l +c [MiUo,Un+l) + MiUn+uUo))"" = -ImK- 

Keeping on the left hand side only what is still unknown, we end up with: 


+ Vm • V^xCTn^l + Cm {(Jm + CTn+l ^8^ 

[ f ;)^ + L{d) (/ - Pm) am Vm) 


Cm 

+M{am Cm, (/ - Pm) U^+x) P M{(1 - Pm) U^+X.^m Cm) + 


where the constant Cm is the one defined in (2.12). We have thus found that the crJJYi’s must satisfy 
decoupled incoming transport equations in x (M/0Z) with infinitely smooth coefficients and source 
terms (all vanishing for t < 0). The only task left is therefore to determine the trace of each m gZ, 

on Xd = 0. 

We recall the decomposition (2.4) introduced in Section 2. We thus introduce a decomposition 


(4.10) 


^ al^^^{t,y,0,eo)rm = an+x{t,y,0o) e+ Un+x{t,y,9o), 

m£X 
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where Qn+i is an unknown scalar function (with zero mean value with respect to 9q), and Un+i takes its 
values in E®(r, r/). Thanks to the compatibility condition (4.6), the function Un+i is uniquely determined 
by solving 


BUn+i - a,-' <9^=00 

VmGX 

- B {{Cide)-^ - (Cide)-^ F:i^)\,^=o,e,=-=e^=9o ) + (G„ -^)(t,y,0o) = 0, 

where dg^ denotes the inverse of dg^ when restricted to zero mean value functions. We get Un+i G 
C'^((-oo°r];^+'^(K'='-^ X (M/0Z))), and W„+i|i<o = 0. 

Before going on, we observe that no matter how we define the scalar function On+i, our definitions so 
far of all components of Un+i give the relation 

(4:-ll) B hln+l\xa=0,ei=---=eN=9o + Gn = 0, 

which means that (4.3) (b) will be satisfied if n = 0 , or (4.3) (c) will be satisfied for fc = n — 1 if n > 1 . At 
this stage, it only remains to determine the scalar function On+i, and the trace of on {xd = 0 } will 
be obtained by computing 

(f, y, 0, ^o) — Un+l (^) yj ^o) T ^n+l,m(^! , 

where Un+i,m stands for the component oiUn+i on in the decomposition ( 2 . 1 ). 

Step 5: determining Part II. Evolution equation for an+i- 

This step mimics the analysis in Section 2, and more specifically Steps 4 and 5 there. Let us assume 
that hln +2 bas no outgoing mode, which will be fully justified at the next order of the induction process. 
The nonpolarized components of the corrector obtained by looking at (4.1) (c) for k = n. For 

m € I, we thus define (/ — Pm) ^^2 ^be solution to 

(/ - Pm) d0^U^_^2 + Rm L((9)(a;(Yl ^rn) + + ^^n+l dj^Pm Rm (dAj(O) • Tm) Tm 

= -Pm K - Rm L{d) (/ - Pm) " M { Un+1 , Um ^m) 

- M{am rm, {I - Pm)U!^+l) - M{{I - PmW^^^am Vm) " (>l(Wo, J , 

where we should keep in mind that is still not fully determined, but the right hand side of the 

equality has already been constructed. Since the partial inverses Pm satisfy Pm Ad{0) Pm = 0, the term 
Pm L(5)(u((Yi ^m) only involves tangential differentiation with respect to the boundary {xd = 0 }, so 
we can take the trace of the latter equation on the boundary. We then substitute a„+i(f, y, 0m) Cm + 
Un+i,m{t,y,0m) for Ud=o ?’m- These operations yield 

(4.12) bB j;(/-Pm)(50^Wr+2)lx,=o.^^^=0o 

m£X 

= -^LopUn+i + (2u - bd^b{0) ■ (e,e)) (aSeoUn+i + Un+i dg^a) + yi,„(f, y, 0o), 

with V and Alop defined in (2.23), and gi^n is explicitly computable from all previously determined 
quantities. 
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We now determine using (4.1)(c) for k = n, computing the noncharacteristic components and 

by taking the trace on = 0. All these operations lead to the equation 


^d8)Kl2\x,=0 = -Yl +-M(^m^^r+1,^0)]“L=0 

m^X 

m^X 

r 1 ric 

- [MiJAoMT+i) + M{KIM\ U.=o - lF“U,=o - mKliUa=o 

We then use the decomposition: 

Pm ^n +1 |3:d=0 — Un+1 (^) Vi 9m) ^m A kin-\-\^m{pi Vi 9m) ) 
where Un+i,m has already been determined, and we thus obtain the expression: 

(4.13) 1x^=0,6i='"=6iv=^'o ~ ®per(U)Un+l) + 92,n j 

where 


Bper[u,u]:=- Y Y1 

7711^7712 k£X 

m\ ,7712 GX 


/ 


\ 


'^km2 A ^fcm2 '^km^ ) L{kmi + km2 dipm2) {km2 Em2,mi) 

^7712 “h^7712 ~ky 
y kmi fc77i2#0 


2i7r kOo/& 

^ 5 


the vectors Emi,m 2 ure defined in (2.25), and g 2 ,n is explicitly computable from all previously determined 
quantities. 

We now consider (4.3) (c) for k = n, which we rewrite equivalently as 


B kin+2 |a:ri=0 4“ B ^ ^ 


77,-1-2 x^d — O,07Ti—^0 


TTl^X 


PB 'Y) Pm)i^n+ 2 \xd= 0 ,dm= 9 o P Bh{j^_^ 2 \xa=O,ei=---= 0 pi= 0 o P '^n+l d 6(0) • (e, c) 

mgl 


(4.14) 


+d^6(0) • ( a e, Un+i + Un+i P ^ ^ {I — Pm) 

\ ttiGX 


+G„ = 0. 


We differentiate the latter equation with respect to 9q, apply the row vector b and use (4.12) and (4.13) 
to derive the governing equation for a^+i: 

(4.15) 2 U (fl S^qO) ALopO-n+l 4“ Qper I®) Un+l] 4“ Qper [Un+l) u] — 9n i 


where g-n incorporates all contributions from the source terms gi,n) 92 ,n and the one obtained after differ¬ 
entiating the last line of (4.14) and applying b. Moreover, Qper is the operator defined in (2.24). 

Observe that the above governing equation (4.15) for a^+i is a linearized version of (2.22). Well- 
posedness for (4.15) follows from the same arguments as those we have used in Section 3, namely from 
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Theorem 3.1 which shows that (4.15) is a transport equation for a^+i that is perturbed by a nonlocal 
’’zero order” term. We thus construct a solution an+i G C°°((—oo, T]; x (M/0Z))) to (4.15) 

that vanishes for t < 0. 

Step 6: conclusion. 

We have now determined each component oilAn+i, which gives a profile in C°°((—oo, T];x 
(M/GZ)-^)) with its spectrum in Zj . Moreover, Un+i vanishes for t < 0, and it satisfies (4.9), (4.11). It 
is also a simple exercise to verify that, if we define Fn+i according to (4.7), our construction of Un+i gives 
the compatibility conditions 


V m € X, irn F^+i = 0 , and Fn+i = 0 , 

which is nothing but (4.5) at the order n + 1. Step 5 above also shows that, with G^+i defined as in (4.8), 
we have the compatibility condition: 

Ud=o,6»^=6»o - -S5eo(('C(^6»)"^-C+i)Ud=o,6»i=-=ejv=e'o) +'^eoG'n+i (T 2/, ^o)) = 0, 

mGX 

which is nothing but (4.6) at the order n+1. We have therefore proved that (H.„+i) holds, which completes 
the induction. 

4.3 Proof of Theorem 1.10 

We now quickly complete the proof of Theorem 1.10. The analysis in Paragraph 4.2 shows that there 
exists a sequence of profiles {Un )n>o in C°°((-oo,r];ll+°°(M^ x (M/(0Z))^)) that satisfies the WKB 
cascade (4.1), (4.3), and Un\t<o = 0 for all n G N. Moreover, each profile Un has its 0-spectrum included 
in Z^. The uniqueness of such a sequence also follows from an induction argument, where we use 
the regularity of each profile to justify all computations that we have made in order to construct the 
UnS (successive differentiations, identification of Fourier coefficients, substitution {9i,... 6 q on 

the boundary etc.). Uniqueness of the sequence {Un)n>o then follows from the uniqueness of smooth 
solutions to all amplitude equations we have had to solve: Equation (2.22) and its linearized version on 
the boundary, Burgers equation (2.12) and its linearized version in the interior. All other operations, such 
as the determination of the non-characteristic components, are of ’’algebraic” type and obviously admit 
a single smooth solution. 

To complete the proof of Theorem 1.10, we thus only need to show that the approximate solutions 
built from the sequence {Un)n>o satisfy the error estimates claimed in Theorem 1.10. We thus consider 

N1+N2 / \ 

n=0 ^ ^ ^ 

Let us first consider the boundary conditions in (1.1). Setting = £Ve{t,y) for simplicity, 

we compute 

Ni+N2+1 n N 1 +N 2 

= E E d”'>(o)'fe.■■■.*>=) 

n=l Ui,...,Un=0 

_A^i+A^2H“2 /*1 

+ (JV, + JV, + 1)1 / . (1.,,d». 
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Collecting powers of e and using the WKB cascade (4.3), we get 




if A^i = A ^2 = 0, and 


) = y (1 - s) d^6(e s Uo) ■ {UoMo) y, 0, 




ds, 


U,=o) -e^G (t, y, 

{Ni+N2+l)'^ N 1 +N 2 +I ^ U \ U \ 

= E E ^ E 

n=A^i+A^2+2 1=1 -^ 

ui,... +-^2 

^A^i+A^2+2 

+ (m + N, + 1)! JJ'- »,) d». 

if A^i + N 2 > 0. 

Each profile l^/n belongs to C°°((— 00 , T]; x (M/(0 Z))^)) and vanishes for t < 0, hence it also 

belongs to L°°((—oo,r] x (M/(0Z))-^;L^(R'^“^)) nL°°((—oo,T] x x (M/(0Z))'^) when restricted 
to the boundary {xd = 0}. We thus have uniform bounds with respect to e of the type 

||d^6(0) • y, 0, (fo/e,..., (^o/e)||Loo(( — cxd,T]; L2(Rd-l)) < C , 

and similarly for the above integral remainders in Taylor’s formula. When Ni + N 2 is positive, we thus 
get 

G{t,y,^o/e)\\L^((-oo,T]-,L^Kd-i)) < Ce^i+^ 2 + 2 ^ 

and we can derive the exact same estimate for the L°° norm of the error at the boundary. 

When differentiating with respect to y, each partial derivative gives rise to a factor 1/e, which yields 

\\H'^T'^’^^’^^\xa=o) — G(t, y, (/?o/e)||Loo((_oo^7^].j^]V2(Rd-i)) < 

The case Ni = N 2 = 0 is similar, except that the error is as large as the source term G, which is not 
so interesting from a practical point of view. 

We leave to the interested reader the verification of the error estimate in the interior domain, which 
involves a little more algebra but no additional technical difficulty. 


4.4 Extension to hyperbolic systems with constant multiplicity 

The extension of the derivation of the leading amplitude equation (2.22) to hyperbolic systems with 
constant multiplicity is not entirely straightforward for two reasons, one related to the zero mean property 
of Uq and the other to the solvability of the interior profile equations. 

When we analyzed the WKB cascade (4.1), (4.3), the first point was to prove that the leading am¬ 
plitude Uq necessarily has zero mean. This property does not extend obviously to the case of hyperbolic 
systems with constant multiplicity (similar issues arise in [CGWll]). Let us recall that when the system 
is hyperbolic with constant multiplicity, Lemma 2.1 and Lemma 2.2 hold. We then use the notation 
introduced in Paragraph 2.1 for the projectors PmiQm and the partial inverses Rm- The only difference 
with the strictly hyperbolic framework is that we do not use the row vectors ^rn here. We now analyze 
the equations (2.6), (2.8) in this more general framework. 
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Equation (2.6) (a) shows that the leading prohle Uq can be decomposed as 

M 

Uoit,X,ei,...,eM)= Uoit, x)+Y, Kit, X, 9m) , Pm K = K, 

?n=l 

where each is ©-periodic and has zero mean with respect to 9m- The mean value Uq satishes the 
homogeneous boundary condition (2.11), and we are going to show that (2.10) still holds. Indeed, the 
mean value Uq satisfies 

L(d)llo + MiUo,Uo) = 0, 

which, in view of the decomposition of Uq and the definition (2.7) of M, can be first simplified into 


M 


muo + y^ M{u^,us^) = o. 


m=l 


We then compute 
(4.16) 

= d^^m (d7l,(0) -U^) deJA^ = d,^m dV,(0) • deJA^) = ^ d,^m de^d^m ■ iK,K) > 

where we have used Assumption 1.2 and the symmetry of d^/j(0). Each Ai{lA^,lA^) therefore has mean 
zero and Wq satisfies both (2.10) and (2.11) as in the strictly hyperbolic case. The result of Step 1 in 
Paragraph 2.2 thus extends to conservative hyperbolic systems with constant multiplicity. 

We derive the interior equation for each IAq^ by retaining only the 0m-oscillations in (2.6) (b) and 
applying the projector Qm- Using [CGIO, Lemma 3.3] and the absence of resonances, we get 


(4.17) 


{dt + v™ • V,) QmK + ^ de^Qm d^fj{0) • iK,K) = 0 ■ 


Let us also recall that the restriction of Qm to Im Pm is injective, so that QmU^ uniquely determines 
lA^^ with lA^ = PmiA^. In spite of the symmetry of d?f{0), (4.17) is not obviously a symmetric hy¬ 
perbolic problem for the unknown Qm^^, so its solvability is not immediately obvious. We encoun¬ 
tered a similar difficulty in [CGWll], and we resolve it here in a similar way using the expression 
djifm (dAj(O) • lA^) dg^lA^ for the nonlinear term and the following lemma. 

Lemma 4.2. Let w G and write d(pm = iz,-q,(Xm) = (-Afc^(0,0, where Xk^{C) = Xk^{u,^)\u=o. 
Then there holds 


(4.18) 


• W I Pm=[- d„Afe^(0,^) • w) QmPm ■ 


i=i 


Proof. For u near 0 let Pmiu) be the projector on 


Ker -Xk^{u,Q I+ yL A 


i=i 



in the obvious decomposition of (that is the analogue of the first direct sum in (2.2) for u close to the 
origin); thus the projector Pm in Lemma 2.2 is Pm = Pm{d). Differentiate the equation 


-Kiu,Ql+ yiyjiy) \ Pmiu) = d, 


i=i 
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with respect to u in the direction w, evaluate at u = 0, and apply Qm on the left to obtain (4.18). □ 

Taking w = in (4.18) and using we see that (4.17) is a symmetric hyperbolic system 

(essentially scalar) for the unknown so it can be solved with appropriate boundary and initial 

conditions just like the corresponding equations in the strictly hyperbolic case. 

Equation (4.17) shows that all outgoing modes in IAq vanish, and there holds 

(4.19) Uo{t, X, 01,..., 0m) = X, 0™), . 

m£X 

In particular, there exists a scalar function a that is 0-periodic with zero mean, such that 

(4.20) 7/,5"(t,y,O,0o) = a(t,y,0o)em, where = T’m e . 

Our goal is to derive the amplitude equation that governs the evolution of a on the boundary. 

At this stage, the analysis of Steps 3, 4, 5 in Paragraph 2.2 applies almost word for word, with obvious 
modifications in order to take into account the possibly many incoming and outgoing phases. Namely, we 
can show that the first corrector Ui has no outgoing mode. We can also determine the non-characteristic 
component and the non-polarized components (/ — P^)1x^=0 on the boundary in terms of 

the amplitude a. We end up with the exact same equation as (2.22), with a real constant v and a vector 
field Xlop as in (2.23). The new expression of the bilinear operator Qper reads (compare with (2.24)): 


(4.21) Qper[a,a] 

( 


m£0 m\<m 2 

mi ,m 2 


E 

fcez 


E 


^mi ~\~^m2 — 
fcmi ^mo 


kmibBAfii^O^ Qm Buii_^iYL2 ^m2 bB A(^( 0 ^ 

kmi (^mi ~ ^m) + ^rn2 (^mo ~ ^m) 




^2in kOo /© 


with vectors as in (2.25). 

The analysis of the WKB cascade (4.1), (4.3) proceeds as before, taking into account that incoming 
amplitudes 7/™ are propagated in the interior domain by vector-valued Burgers type equations (4.17) and 
appropriate linearizations at which can be solved using Lemma 4.2. 


36 





Part II 

Pulses 


5 Construction of approximate solutions 


We follow the approach of Section 2 and first deal with strictly hyperbolic systems of three equations. 
We keep the notation of Paragraph 2.2, and make Assumption 2.3. Let us now derive the WKB cascade 
for pulse-like solutions to (1.1). The solution Ug to (1.1) is assumed to have an asymptotic expansion of 
the form 


(5.1) 


Ue ^ e'^e^Uk 
k>0 



£ 



We use the notation 9q as a placeholder for the fast variable <y9o/e, and for Xd/£- Plugging the ansatz 

(5.1) in (1.1) and identifying powers of e, we obtain the following first three relations for the lAkS (observe 
the slight differences with Paragraph 2.2, though we keep the same notation): 

(a) /:(50o,%)ZYo = 0, 

(5.2) (b) C,{de,,dQUi+L{d)Uo + M{UoMo) = ^, 

(c) C{de,,d^,)U 2 + L{d)Ui+M{UoMi)+M{UiMo) + = 0, 

where the differential operators £, A4,A/i,A/2 are now defined by: 

-^(900,%) := Arf(O) (% +iA{T,i^ de^) , 

M(v,w} := djipo (dAj(O) • v) de^w + (dAd(O) • v)d^^w, 

N'i{v,w) := (dAj(O) • v)djW , 

Af2{v,v,w) := (d^^j(O) • {v,v))dg^,w + ^ (d'^AdiO) ■ iv,v))d^^w. 

The equations (5.2) in the domain (—oo, T] xM^ are supplemented with the boundary conditions 

obtained by plugging (5.1) in the boundary conditions of (1.1), which yields (recall B = d6(0)): 

B IAq = 0 , 

BWi + i dH{0) • {UoM) = G{t, y,eo), 

where functions on the left hand side of (5.4) are evaluated at Xd = Cd = 0. In order to get Ue|t<o = 0, as 
required in (1.1), we also look for profiles Uk that vanish for f < 0. 

The construction of prohles in the wavetrain setting was accomplished by first assuming that solutions 
exist within the class of periodic functions of 6 , and then using that assumption to actually construct 
periodic solutions. The construction of profiles in the pulse setting has the difficulty that it is not so 
clear at first in what function space(s) solutions should be sought. Construction of each successive pulse 
corrector involves an additional integration over a noncompact set. Thus, each corrector Uj is “worse” 
than the previous one and successive correctors must be sought in successively larger spaces. In 

fact we will find that correctors beyond U 2 are useless; they grow at least linearly in {6o,^d), and are thus 
too large to be considered correctors. 


(5.4) 


(a) 

(b) 
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We now define spaces Vp C Vh C C^, where is the space of functions K(t,x,6Q,^d)^ valued in 
and bounded with their first-order derivatives. The variables {t,x) lie in (—oo,r] x while the 
variables He in M x M+. In the subsequent discussion we will assume Uq € Vp, Ui € Vh-, and 

U 2 € cl are solutions to (5.2), (5.4), and then construct prohles with those properties. 

Definition 5.1. (a) Let Vp denote the space of functions 

3 

F{t,x,eo,^d) = '^Fi{t,x,9o,Cd) Ad{ 0 ) n , 

i=l 

where eaeh funetion Fi is a finite sum of real-valued functions of the form 

f{t,x,eo+^k^d), f{t,x,eo + ^Q)g{t,x,9o +uim^d), 

and f{t, X, 9o + U) g{t, x, 9o + ^ U) x, 9o + ^d) , 

where the indices k, r lie in {1,2,3}.® The funetions f{t,x,9), g{t,x,9), ete., in (5.5) are and 
decay with their first-order partials at the rate 0{{9)~‘^) uniformly with respect to {t,x). We refer to these 
funetions as the “eonstituent funetions” of F £ Vp. 

(b) Define a transversal interaction integral to be a function ^(t, x, 0O) ^d) of the form 

r+°° 

(5-6) limits X, 9o, U) = - or{t, x, 6»o +kliU + kli) s) T{t, x, 9o +kli^d + [klm - kli) s) ds , 

->id 

where uh, and Wj are mutually distinct. The functions a{t,x,9), T{t,x,9) are real-valued, and 
decay with their first-order partials at the rate 0{{9)~^) uniformly with respect to (t,x). 

(c) Let Vh denote the spaee of funetions 

3 

H{t,x,9o,^d) = '^Hi{t,x,9o,U)Ad{0)ri , 

where each Hi is the sum of an element of Vp plus a finite sum of terms of the form 
(5.7) Il^{t,x,9o,f,d) or a{t,x,9o + klkU) Jp,q{t,x,9o,f,d) , 

where and Jp^ are transversal interaction integrals and the indiees i,l,...,q lie in {1,2,3}. The 
function a{t,x,9) is real-valued, and decays with its first-order partials at the rate 0{{9)~^) uniformly 
with respect to {t,x). 

(d) For H € Vh we can write H = Hp + Hj, where Hp £ Vp and Hj ^ Vp has components in the 
span of terms of the form (5.7). The “constituent funetions” of H include those of Hp as well as the 
functions like a, a, r as in (5.6), (5.7) which constitute Hj. 

Remark 5.2. The same spaees of functions Vp and Vh are obtained if one begins by writing 

3 _ 3 _ 

F{t,x,9o,^d) = '^Fiit,x,9o,^d)ri , H{t,x,9o,^d) = '^Hi{t,x,9o,Cd)ri , 

2=1 2=1 

and imposes the conditions in (a) (resp. (c)) on the Fi (resp. Hi). 

^Elements of Vf with terms involving no triple products were called functions of type F in [CW13, CW14]. 
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5.1 Averaging and solution operators 

To construct the profiles lAQ,Ui,hl 2 , we must solve equations of the form C{dgQ,d^^)U = H, see (5.2), 
where H lies in Vh and sometimes in Vp- In this section we define averaging operators Ep and Eq and a 
solution operator i?oo (involving integration on a noncompact set) that provide a systematic way to study 
such equations. Henceforth we write C := C{dgQ,d^^). The following simple lemma implies the existence 
of most of the limits and integrals that appear below. 

Lemma 5.3. Let a{t,x,6), T{t,x,6) be continuous functions such that 

\a{t,x,6)\ + |r(t,x,0)| < C {6)~^ , 


for some N > 2, uniformly with respect to {t,x). Let i,i,m,q lie in {1,2,3} and suppose i,£,m are 
mutually distinct. Then the following estimates hold true uniformly with respect to all parameters: 



|cr(t,X,(9o + + (W£ 

a{t, X, 00 + Wj Cd + (w^ 


wJ'S)|ds < 1, 

Wj S) r(t, X, 00 +^ifd+ (Wm 


Wi)s)|ds < (^d) 


r-\-oo r-\-oo 

{c)forN>3, / \cr{t,x,0o+^fd + i^i 

hd Js 




T{t, X, 00 +^gfd + i^i - Wg) s + (Wm - ^i) r) I dr ds < {(d) ■ 


Proof. Part (a) is immediate. To prove (b), we use the Peetre inequality: 


r+oo 

/ |cr(t, X, 00 +^i(d + (^£ - Wj) s) T(t, X, 00 +^i(d+ (Wm “ Wj) s)! ds 
hd 

P + OO /•+00 

^ / (6'o + Wi?d + (w£-wJs)“^(6'o + WiCd + (wm-Wi)s)“^ds < / ((w^-a;^)s)“^ds, 

hd hd 

and the result follows (recall > 0). Part (c) is also proved by using the Peetre inequality: 


p+oo P+OO 

/ / W{t,x,0o + Ug(d + {^i-lilq)s + {ui^-iwf)r)- 

J^d ds 

p-\-oo p-\-oo 

r{t,x,0o + i^(d + i^i-^q)s + {uim-^i)r)\drds< / ((^^ - dr ds < 

J£d J s 


-N+2 


□ 


The definition of Vh and Lemma 5.3 imply that the limits and integrals in the next definition are well- 
defined. 


Definition 5.4 {Ep, Eq, Roo)- For H G Vh, define averaging operators 

EQH{t, X, 00 , (d) := f J ^0 + {(d - s), s) ds^ Hrf(O) Vj , 

EpH{t,x,0o,(d) J ij ■ Ad{O)H{t,x,0o+uij{(d- s),s)ds'^ rj . 
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For H E Vh such that Eq H = {), define the solution operator 


Roo E ^5 ^0 ; ^d) 



Ij ■ H{t,x,eo + Uj {Cd 


s), s) ds I r 


Remark 5.5. (a) Suppose F S Vp reads 


3 

F{t,x,eo,^d) = '^Fi{t,x,9o,^d) Ad{0) ri, 

i=l 


where each Fi has the form 

3 3 

Fi (i, X, ^0? f,d') ^ ^ ri}, ff^(t^ X, ^0 F f,d) “1“ ^ ^ ^i,m ^0 F 0^ ^d^ (^? ^0 d“ 5 

k = l 1^771=1 

Then 

3 

EqF = Y1 (“i ^0 + F-i U) F bF gl(t, X, 9o + U) Kit, X, 9o + Wi (d)) Ad{0) ri. 

i=l 

The obvious analogues hold when F G Vp involves triple products, or when Ep is used in place of Eq. 

(b) Let H G Vh and write H = Hp + Hi, where Hp G Vp and Hj ^ Vp as in Definition 5.1(d). Then 
Lemma 5.3 implies Eq Hi = Ep Hi = 0. 

The following Proposition gives the main properties of the operators Ep, Eq, Roq. 

Proposition 5.6. Suppose H G Vh and recall the notation L = £{ 80 ^, 8 ^^)- Then there holds: 

(a) EqCH = CEpH = 0. 

(b) If Eq H = 0, then Roc H is bounded and C Roc H = H = {I — Eq) H. 

(c) IfEpH = 0 ,then Roo £ H = H = {I - Ep) H. 

Proof, (a) £EpH = 0 follows directly from Remark 5.5 (a),(b). To show Eq £ H = 0, we write H = 
Hi ri, and note that 

£ {Hi n) = (% - Wj 000) Hi ^d(O) ri. 

Thus, the integrals in the definition of Eq £ H can be evaluated and are uniformly bounded with 
respect to T. After dividing by T, the limit as T goes to infinity vanishes. 

(b) Boundedness of RooH follows from Lemma 5.3. Writing H = Hi Ad(0) rj, a direct computa¬ 
tion using L{dipi) = 0 shows £ Roo H = H. 

(c) With /? = (r, 77 ), define A{f3) by L{dipk) = A{/3) Fui^AdfA), and observe that 

(5.8) 4 A(/3) rj = -Wfc djk . 

Writing H = Hj rj and using (5.8), a direct computation yields 

^ ^ / r +00 ^ \ _ 

Roc£H = - (% - Uj 80 ^) Hj{t, X, 9q + uij (4 - s), s) ds ] 4 = X] ^1 ’ 

i=i ^ j 

since limjvf^.-i-cxD Hj{t, x, 9 q + Uj (4 — M),M) = 0 when Ep H = 0. □ 
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The next proposition, which extends [CW13, Proposition 1.22], is used repeatedly in constructing 
IA\ and IA 2 hereafter. 

Proposition 5.7. Suppose H G Vh- 

(a) The equation CU = H has a bounded solution if and only if Eg H = 0.® 

(b) When Eg H = 0, every bounded solution to the equation CU = H has the form 

3 

U = ^ X, 6q T T RooH. 

m=l 

(c) When H € Vp satisfies Eg H = 0, solutions to CU = H satisfy 

3 

EpU = ^ Tm{t,x,eQ + u^id)rm , {I “ Ep)U = RooH . 

m=l 

Proof, (a) The direction -4= is given by Proposition 5.6(b). (=>) Suppose there is a bounded solution to 
CU = H, and write H = Hp + Hj as in Remark 5.5(b). Since Eg Hj = 0, there is a bounded solution 
io CU = Hi, so we conclude there is a bounded solution to CU = Hp, and similarly io CU = EgHp. 
From the explicit form of Eg Hp given in Remark 5.5(a), we see that there is no bounded solution of 
CU = EgHp if EgHp 7 ^ 0 (for otherwise solutions display a linear growth with respect to ^rf). Thus 
Eg Hp = 0 and hence Eg H = 0. 

(b) By Proposition 5.6(b), RooH is a bounded solution oi CU = H; moreover, the general 

solution oi CU = 0 has the form X]m=i '^rn{t, x, 9 q + I’m- 

(c) Lemma 5.3 implies Ep Roo H = 0 when H G Vp and Eg H = 0. □ 

5.2 Profile construction and proof of Theorem 1.11 

For Qp := (—oo,T] x x we define the weighted Sobolev spaces: 

r^(Pr) := {u G ^^(Pr) : 0^ G ^^(Pt) if a + |/3| < A:} . 

This is a Hilbert space for the norm 

a+\P\<k 

and r^(P 7 ’) is an algebra for k > {d + 2)/2. For L G N, when k > + L + 1, elements of r^{Q,T) are 

and decay with their first-order partials at the rate {0)~^ uniformly with respect to {t,x). 

Definition 5.8. (a) Suppose k > + 3. Let Vp C Vp be the subspace consisting of elements whose 

constituent functions lie in r^(P'r). 

(b) Suppose k > -|- 4. Let V|) C Vh be the subspace consisting of elements whose constituent 

functions lie in 

®The proof actually shows that CU = H has a W solution that is sublinear with respect to if and only if Eq H = 0. 
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For Kq > 8 + we assume now that the equations (5.2), (5.4) have solutions Uq G Vi?, Ui G Vh-, 
U 2 G cl such that 


(5.9) 

(a) Uo G , Uo 


Ep IAq 


3 

^ am{t,x,eQ + u^£,d)rm , 

m=l 


where £^(t,x) := / fTm(t, a:, 0^) d6»m = 0 , 

JR 


(6) EpUi G and (/ - Ep)Ui G , 


and then construct solutions with those properties. 


Step 1; the leading profile Uq has no outgoing component. 

This step justihes one of the causality arguments used in [MR83]. Equation 5.2(a) and Proposition 
5.7(b) imply that the expression of Uq reduces to: 

3 

UQ{t,x,0Q,id) = ^ am{t,x,eQ + u^id)rm 

m=l 

for functions Um to be determined. The last variable of am is denoted 9m in what follows. 

Since Ui is a bounded solution to CUi = where Fq = —L{d)UQ — AA.{Uq^Uq) G Proposition 

5.7(a) implies EqFq = 0; that is, 

(5.10) + + = m = l,2,3, ^ ■ r„)r^ ^ 

Since ip 2 is outgoing, this implies <72 = 0, and the boundary condition (5.4) (a) gives, as in Paragraph 2.2, 
the existence of a scalar function a such that 

(5.11) cri(t,y,0,6»o)ri = a(t,y,6»o)ei , cr3(t, y, 0,6*0) 73 = a(t, y, 6*0) 63 . 


Step 2: showing (/ — Ep)Ui G 

At this stage, we know that the leading profile Uq reads 

(5.12) Uo{t,x, 6 o,U) = cri{t,x,eo + ^d) n + a3{t,x,0o + , 

where ai,a 3 satisfy (5.10) and their traces satisfy (5.11). We thus compute 

(5.13) Fo = -L{d) (cJi ri + CJ 3 73 ) - fJi dg^ai i?iq - 0-3 d 0 ^a 3 ^ 3,3 - 03 de^ai Ri ^3 - ai 80^03 R 34 , 
with 

V TTii, ?7l2 — 1)3, Rmi ,m 2 • — djfmi (dA^ (0) • Tm 2 ) Xmi ) 

and the functions (cJi, de-^ai), resp. (cr 3 , do^a 3 ), in (5.13) are evaluated at (t, x, 0o + iiii Cd), resp. (t, x, Oq + 
W 3 ^rf). By Proposition 5.7(b), we have 

3 

^ ^ Xm{t, X, 9 q F F.m ^101 T Eoo Fq , 

m=l 

3 / ^+oo 

(/ - Ep)Ui = RooFo = - ( / Fm{t,x,9o + {^d - s),s) ds 

m=l 
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Since we have Eq Fq = 0, the latter integrand x, 9 q + — s), s) reads 

^k{t, X, 6*0 + + (wfc — ^m) 5 )+ 

k^m 

ak{t, X, 6*0 + + (Wfc - s) 90 ^crfc(t, X, 60 +^m^d + i^k “ ^m) s) + 


k^m 


(5.14) 


Y ^0 + -"1 de^(^e{t, X, 00 + Wm Cd + - Wm) s) + 

l^m 

Y ^0 + iiim Cd + (lii£ - yim) s) dg^Umit, X, Oq + ^rf) + 


i^m 


Y X, 00 + Wm Cd + (Wf - <^m) s) dg^akit, X, 0o + Cd + (Wfc “ ^m) '«) ; 

i^k^i^m^k^m 

where the c™, are real constants, and is the tangential vector field given by 

14 ™ (Tfc := im L{d) {akVk) ■ 

After integrating with respect to s on [^d, + 00 [ the second and third lines of (5.14), constitnent fnnctions 
in Vp° are obtained. The first and fourth lines can be integrated using the moment zero property of the 
cjm’s to yield constituent functions in Vp°~^ and Vp°~^ respectively.^*^ The integral of the fifth line is a 
linear combination of transversal interaction integrals with constituent functions in 

Step 3: equations for EpUi. 


Since U 2 € is a solution to CU 2 = Fi, where 


(5.15) 


:=- 


L{d)Ui + M{UM) + M{UiMo) + J^iiUoMo) + muMMo) 


G V 


Ko-A 
H J 


Proposition 5.7(a) implies EqFi = 0. We rewrite this and include the boundary condition on lAi to 
obtain 


(a) Eq L{d) EpUi + M{Uo,EpUi) + M{EpUiMo) 


(5.16) 


-E, 


Q 


L{d) (I - Ep)Ui + M{Uo, (/ - Ep)Ui) + M{{I- Ep)UiMo) +M 1 FM 2 


( 6 ) BEpUi = G - -d26(0) • {UoM^) -B{I-Ep)Ui, onxd = ^d = 0. 


Letting V := (/ — Ep)Ui (which has been constructed from Uq in Step 2 above), and decomposing 
V = Vp + Vi as in Definition 5.1(d), we can use Lemma 5.3(c) to see that 


E, 


Q 


L{d) Vi + M{Uq, Vi) + M{ViMo) 


= 0 . 


Thus, (5.16) (a) simplifies to 

(5.17) Eq \L{d)EpUi + M{UQ,EpUi)+M{EpUiMa) 


= -E, 


Q 


L{d) Vp + M{Uo, Vp) + M{Vp,Uo) + Wi + A 62 


G V 


Ko-3 


10 


Without this moment zero property, these integrals and thus Ui would be no better than bounded; consequently, U 2 
would be unbounded. 
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The components of EpUi will be determined from equations (5.17) and (5.16)(b). 

Step 4: determining the component of EpUi on r 2 . 

We now show T 2 = 0, thereby justifying another causality argument in [MR83]. Since Uq is ’’purely 
incoming” (no {Oq + ^2 ^d)-dependence), using Remark 5.5(a), we see that the component on r 2 of the 
right side of (5.17) is zero. This is rather clear for the terms Eq (A/i + A/ 2 ), because A/i + A /2 is a linear 
combination of terms of the form 

fiit,x,9o + , f3it,x,eo EuisCd), g{t,x,eo+^^^d)h{t,x,eo+U3^d), 

and it is also true for the terms arising in Eq [L{d) Vp + Vp) + M.{Vp,V(o)] after examining the 

form of Vp = {Roo T'o)f- 

Similarly, the component on r 2 of Eq [M.{Uo, EpUi) + Ai{EpUi,Uo)] is zero, whatever the value of 
the constituent functions Tm in EplAi (which remains to be determined). This forces the function T 2 in 
the decomposition of EpUi to satisfy the homogeneous transport equation 

{dt + V2 • Vx)t2 = 0, 

Since (p 2 is an outgoing phase, T 2 is identically zero. 

Remark 5.9. The justification of T 2 = 0 relies on the assumption that {b.2)(c) admits a bounded solution 
U 2 - Boundedness ofU 2 makes the expression e^U 2 meaningful as a corrector to the approximate solution 
eU(j +e^Ui. However, we shall see in Appendix B that assuming boundedness ofU 2 has a major conse¬ 
quence on the leading order profile Uq. In particular, the governing equation (5.19) below for the evolution 
of a on the boundary will not coincide with the equation obtained by considering (2.19) in Part I in the 
regime 0 + 00 . 

Step 5: the Mach stem equation for a. 

From the previous step of the analysis, the trace of the first corrector Ui satisfies 

^+00 

Ui{t,y,0,eo,0) =i^ri - / £2 J'o{t,y,0,do - UI 2 X, X) dX r 2 + ^r^ , 

Jo 

where * denotes a coefficient whose expression is not useful for what follows, and Eq is given by (5.13). 
Plugging the latter expression in (5.4) (b) and applying the row vector b, we get 

^+00 1 

(5.18) -bBr 2 £2Eo{t,y,0,9 o-ui2X,X)dX+ -bd^b{0)-{e,e)a^ = bG, 

Jo 2 

which is the solvability condition for EpLIi\x,i=^,^=o in (5.16). R remains to differentiate (5.18) with 
respect to 9q and to identify the first term on the left hand side of (5.18). More precisely, we compute 

J 2 Eo{t, y, 0, 9o,^d) =-£2 Ttan('9) {a{t, y, 9 q + id) ei + a{t, y, 9o + Wg id) £ 3 ) 

-^£ 2 Eiq dg^{a^){t, y, 9o + Wg id) - ^ ^2 -E/ 3,3 dg^{a‘^){t, y, 9o + Wg id) 

- £2 -E/ 1,3 {dgoa){t, y, 6»o + Wg id) a{t, y, 9 q + Wg id) 

- £2 -E/gq a{t, y, 9 q + Wg id) {dg^a){t, y, 9o + Wg id), 
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with Ei^i, as in (2.25), and 


d-l 

Lts,n{d) ■=dt + '^ Aj{0) dj . 
i=i 

Using the expression of the matrices in Paragraph 2.1, we thus find that (5.18) reduces to 

(5.19) vdooia^) - Xi,opa +deQQpn\[a,a] = bde^G, 

with V as in (2.16), Xuop as in (2.17), and^^ 


^+oo 

(5.20) Qpui[a,a]{9o) := {bBr2)hEi^3 / (9eoa(6'o + (wi - W2) X) a( 6 »o + (^3 - ^2) X) dX 

Jo 

r+00 

+ {bB r2) £2 / o,{9q + — ^ 2 ) E — ^ 2 ) 

Jo 


dX. 


Step 6: completing the construction of Up, Up, U 2 - 

It is proved in Corollary 6.5 of Section 6 that, with Kp, Kp € N and G as in Theorem 1.11, there exists 
T > 0 and a unique solution 

a G nfioC'((-oo,r];r^i-i-^(Rj^,)), 

to (5.19), (5.20). From (5.11), this determines the boundary data of ap, ap. Corollary 6.6 of section 6 
yields fTi,cJ 3 € r^°(U 7 ’) satisfying 5.10 (up to restricting the time T > 0). Moreover, a and thus iTi,iT 3 
are shown there to have moment zero. That completes the construction of with the properties in (5.9). 

Assuming these results for now, we complete the construction of Up and U 2 - With Up determined, 
(/ — Ep)Up G is now constructed as in Step 2. To determine EpUp we return to (5.17) and 

(5.16)(b), noting that the right side of (5.17) is now determined. Writing {I — Ep)Up = Vp + Vj as before, 
we have^^ Vf\x^=^^=o G and the same holds for the trace of V/ as a consequence of Corollary 

6.3. The right hand side of (5.16) (b) satisfies the required solvability condition, so (5.16) (b) uniquely 
determines the trace EpUp\x^=^j=p G taking its value in E®(r, r/); recall (2.4). Equations (5.17) 

and (5.16) (b) determine decoupled transport equations for the components Tp, Tp of EpUp, and we obtain 
G and hence EpUp e 

The profile U 2 satisfies CU 2 = Ep, where Fp as in (5.15) satisfies Eq Fp = 0. Thus, Proposition 5.7 
yields a solution U 2 = Roo Fp G . 

Apart from the results proved in section 6 that were used in this step, this completes the proof of 
Theorem 1.11 in the 3x3 strictly hyperbolic case. The profile Up satisfies (5.2)(a), (5.4)(a), and the 
correctors Up, U 2 satisfy (5.2) (b), (c) and (5.4) (b). 

6 Analysis of the amplitude equation 

6.1 Preliminary reductions 

Our goal in this section is to prove a well-posedness result for the “Mach stem equation” (5.19). We 
focus on the case of 3 x 3 strictly hyperbolic systems, and leave the minor modifications for the extension 

^^The variables {t,y) enter as parameters in the definition of Qpui so we omit them. 

^^We use self-explanatory notation here. 
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to iV X systems to the interested reader (see Paragraph 6.5 for the derivation of the corresponding 
amplitude equation). Up to dividing by nonzero constants, and using the shorter notation 0 instead of 
6 q, Equation (5.19) takes the form 

( 6 . 1 ) Sjo + w • Vya + cadga + deQ^vM-, o\ = g , 

with -w € c G M, and the quadratic operator Qpui is defined by 


r+oo 

Qpui[a,a](6») :=/xi / dga{9 + - U 2 ) X)a{e + - U 2 ) X) dX 

Jo 

r+oo 

+ g 3 a{9 + - U 2 ) X) dga{9 + {(^ - U 2 ) X) dX , 

Jo 


where ^i ,//3 G M and the w^’s are pairwise distinct. The latter operator only acts on the 0-variable, and 
(t, y) only enter as parameters, which we do not write for simplicity. 

We first reduce the expression of Qpui by recalling that , ^3 are the two incoming modes while W 2 is 
the outgoing mode. There is no loss of generality in assuming Then we define the two nonzero 

parameters 

hfs “ hfi his “ hli 

that satisfy (53 = 1 -|- (5i. Changing variables in the expression of Qpui and redefining the constants /ii^ 3 , 
we obtain 


r-\-oo r-\-oo 

Qpui[o,h](0) = //i / dea{9 + 5iX)a{9 + 53 X)dX + yi 3 / a{9 + 5i X) dea{9 + 83 X) dX . 

Jo Jo 


For later use, we define the following bilinear operator Fp^i acting on functions that depend on the 
variable 0 G M (whenever the formula below makes sense): 


( 6 . 2 ) 


]Fpul(^^,'f^)(6') 



+ (5i X)v{9 + 83 X)dX. 


The operator Fpuj satisfies the properties: 


(6.3) 

(6.4) 


(Differentiation) 
(Integration by parts) 


90(Fpui(u,u)) = Fpui(50U,u) -hFpui(u,50u), 




1 

— —uv 
83 


-rTpui((90n,u). 
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Using the properties (6.3), (6.4), we can rewrite Equation (6.1) as 


(6.5) dta + 'w ■ Vya + cadga + g¥pui{dga,dga) = g , 

with suitable constants that are denoted c and /i for simplicity and whose exact expression is useless. Our 
goal is to solve Equation (6.5) by a standard fixed point argument. The main ingredient in the proof is 
to show that the nonlinear term Fpui( 50 a, 5ea) acts as a semilinear term in a suitable scale of Sobolev 
regularity. 
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6.2 Boundedness of the bilinear operator Fpui 

The operator Fpuj is not symmetric but changing the roles of (5i and 63 , the roles of the first and second 
argument of Fpui in the estimates below can be exchanged. This will be used in several places. 

Let us first recall the definition of weighted Sobolev spaces: 

:= |u G x M^) : 6 ^ G if a + |^| < A:} . 

This is a Hilbert space for the norm 

ll^llr''(M‘*) ^ 11^ • 

a+\0\<k 

Following the same integration by parts arguments as in [CW13, Proposition 3.3], there holds 
Lemma 6.1. For all integer k, the space coincides with 

|m G X Me) : G l2(M'^)| , 

and the norm o/r^(M'^) is equivalent to the norm 

11^^ “IIl 2(R'') + lkllH'=(Rrf) • 

Our main boundedness result for the operator Fpui reads as follows. 

Proposition 6.2. Let ko denote the smallest integer satisfying fco > (d + l)/2. Then for all k >2kQ + l, 
there exists a constant Ck satisfying 

(6.6) Vtt, u G r^(M'^), ||Fpui(deM, i9eu)||pfc(jjd) < ||u||rfc(iRd). 

The estimate (6.6) of Proposition 6.2 is not tame, but it will be sufficient for our purpose since we 
shall only construct finitely many terms in the WKB expansion of the solution to (1.1) (opposite to 
what we did in Part I where we constructed approximate solutions of arbitrarily high order). 

Proof. Let us first observe that when u,v belong to the Schwartz space 5(M'^), Fpui(u, u) also belongs 
to 5(M'^). By a density/continuity argument, we are thus reduced to proving the estimate (6.6) for 
u,v € S{W^). The decay and regularity of u,v will justify all the manipulations below. 

1) We start with the basic estimate of the function Fpui(M, u). Using Cauchy-Schwartz inequality, 
we have 

+00 2 

u{y, 9 + 61 X) v{y, 6 + d^X) dX 

r+oo /*+oo 

< / \u{y,e + 6 iX)\dX / \u{y,9 + 6 iX)\\viy ,6 + 63 X)\‘^dX 

Jo Jo 

<T^ [ \u{y,e')\de' [ \uiy,e + 6 iX)\\viy,e + 63 X)\^dX. 

I<^i| Jr Jr 
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Integrating with respect to {y, 6 ), and changing variables (use ^3 — di = 1), we get 


l|lFpui(M,u)||^ 2 (Rd) < \^JJu{y,e)\dej \^JJv{y,e)\ d 0 J dy 


< 


TT 


il + 9^)\uiy,e)\^d9] / \v{y,9)\^d9] dy 


l<^l| 7 r‘^-i 

- 1^ ( / (1+ I ||u||^2(Krf). 


Since ko — 1 > {d — l)/ 2 , we have 


i{y,9)\^<C ^ [ \d^u{y,9)\^dy, 


|q:|<A:o —1 ' 

by Sobolev’s inequality, and we thus get (with a possibly larger constant C) 


(6-7) ||Fpui(u, u)||^2(Rd) < d ||u||pfeo(]gd) ||u||2,2(Rd) . 

The ’’symmetric” inequality 

(6-8) l|IPpul(^^5 ^^)||L2(Rd) < C ||^i||l,2(Rd) ||u||pi;o(Rd) , 

is obtained by exchanging the roles of 5i and 83 as explained earlier. 

2) Let us now estimate the Ll^-norm of ¥py,\{dQU, dgv) with k > 2 ko + 1. We first apply the estimate 
(6.7) for the L^-norm: 

||lFpui(96i?^) dgv)\\] 32 ^^d'j < C ||tt||r''o+i(M‘*) ^ ll^llr'=(R'i) ll^llr''(K'i) • 

Using Plancherel’s Theorem, it is sufficient to estimate the fc-th derivatives of Fpui(clett, in order to 
estimate all derivatives of order less than k. Let us therefore consider a multiinteger a of length k, and 
apply the Leibniz formula (this is justified because the differentiation formula (6.3) holds not only for the 
0 -derivative but also for the y-derivatives): 

(6.9) (9" Fpui(50U, dgv) = ^eu, dgv ), 

/3<q: 

where * denotes harmless numerical coefficients, and stand for possibly mixed y,9 derivatives. 

We begin with the extreme terms in (6.9). If |/3| = 0, we need to estimate the term ¥p^i{dgu,d°‘ dgv) 
which, using (6.4), we write as 

-^dgud'^v - ^Fpui(9|0U,5“u). 

We get the estimate 

||Fpui(90U,9"50u)||p2(Rd) <c\\dgu\\ L°°(K'i) ll^llH'=(Rd) + C* ||tt|lr'' 0 + 2 (Md) ll'f^ll//fe(R<i) ) 

where we have used (6.7). Since k > d/2 -|- 1 (this follows from the assumption k > 2A:o -|- 1), we can 
apply Sobolev’s inequality and get 

||Fpul(96lU, 5“ 96l^^)|lL2(Rd) < C ||u||pfc(Rd) ||u||pfc(Rd) . 
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The second extreme term Fpui(5“ Squ, Sqv) is dealt with in the same way. 

Using the assumption A: > 2 /cq + 1, we verify that for (3 < and /3 7 ^ 0, /3 7 ^ a, one of the following 
two properties is satisfied 

(|/3| > 1 and 1 + 1/3| < A: — ko) or (|a| — 1/3| > 1 and |/3| > A;o + 1). 

In the first case, we use (6.7) and get 

||lFpui(9^90U,(9“"^50n)||i2(Kd) < C \\d^ dgu\\-pko(^^d) dev\\L2(^^d-^ < C ||n||r;=(Kd) H'f^HrUR'*) ’ 

and the second case is dealt with in a symmetric way (using ( 6 . 8 ) rather than (6.7)). 

3) It remains to estimate the L^-norm of 9 ^¥p^ji{dgu,d0v). We write (use ds — di = 1 again) 

0 = da (0 + di s) — hi {6 + ha s), 


which gives 


k 

( 6 . 10 ) 9^ ¥p^i{d0U, d0v) = ^ dgu, 9’"~^ d0v ), 

1=0 

with, again, harmless binomial coefficients that are denoted by *. Let us first consider the extreme terms 
in the latter sum and focus on Fpui(cl 0 U, 9^ dgv). We write 

9’^d0v = d0{9^v)-k9^-^v, 

and use the property (6.4) to get 


Vpvii{d0U,9^ dgv) = -^dgu{9’^v) - ^¥p^i{dggu,9’^v) - k¥p^i{dgu,9^ S). 

The L^-estimate follows from (6.7) and from the Sobolev imbedding Theorem: 

||lFpui(i96)U, 9 ^ 90u)||i^2(]Rd) < C ||tt||pfc(Rd) ||u||pfc(-Rd) . 

The estimate for Fpui(6*^ dgu, dgv) is similar. 

For 1 < j < A: — 1, we observe again that there holds either j < A: — A:o — 1 or j > A;o + 1, so that we 
can directly estimate all intermediate terms in the sum (6.10) by applying either (6.7) or (6.8). The proof 
of Proposition 6.2 is now complete. □ 

As an immediate corollary of the proof we have: 

Corollary 6.3. Let kp denote the smallest integer satisfying kp > {d + l)/2. Then for all k >2 kp, there 
exists a constant Ck satisfying 

yu,v € r^(M'^), ||Fpui(u,u)llpfc(Kd) < Ck ||'«||r'=(E'i) lbllr'=(R‘*) • 
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6.3 The iteration scheme 

In view of the boundedness property proved in Proposition 6.2, Equation (6.5) is a semilinear perturbation 
of the Burgers equation (the transport term w • Vy is harmless), and it is absolutely not surprising that 
we can solve (6.5) by using the standard energy method on a fixed point iteration. This well-posedness 
result can be summarized in the following Theorem. 

Theorem 6.4. Let ko be defined as in Proposition 6.2, and let k > 2kQ + 1. Then for all R > 0, there 
exists a time T{k,R) > 0 such that for all data oq € r^(M'^) satisfying ||ao||r'“(R<*) ^ there exists a 
unique solution a G C([0, T]; r^(M'^)) to the Cauchy problem: 

{ dta + 'w ■ VyU + cadea + /xFpui(56ia, 5ea) = 0, 
a|f=o = «o • 

Of course, one can also incorporate a nonzero source term g G T^([0, Tq]; r^(M'^)) in the equation, and 
obtain a unique solution o G C([0, T]; r^(M'^)) on a time interval [0,T] with T < Tq. We omit the details 
that are completely standard. 

Proof. We follow the standard energy method for quasilinear symmetric systems, see for instance [BGS07, 
chapter 10], and solve the Cauchy problem by the iteration scheme 

Idta’^+^ + w • Vya^+^ + ca^ dsa^+^ + fiFp,,i(dsa^,d 0 a^) = 0 , 

|t=0 ®0,n+l ) 

where (ao,n) is a sequence of, say, Schwartz functions that converges towards oq in r^(M'^), and the scheme 
is initialized with the choice = ao,o- Given the radius R for the ball in r^(M'^), we can choose some 
time T > 0, that only depends on R and k, such that the sequence (a”) is bounded in C([0,Tj;r^(M'^)). 
The uniform bound in C([0, Tj; 77^(M'^)) is proved by following the exact same ingredients as in the case 
of the Burgers equation, and the weighted bound is proved by computing 

dt(6^ a”+^) + w • Vy(0^ a"+^) Tea” 5^(0^ a^+^) = -/r 6^ Fpui(50a”, dea^) + kca^ {9^-^ dea ^+^), 

and performing the usual L^-estimate for the transport equation. 

It remains to show that the iteration contracts in the C([0, Tj; L^(M'^))-topology for T small enough. 
This is proved by defining := and computing 

+ w • + ca" = -cr” - //Fp,l(a,r^ 9,a") - /iFp,i(a,a"-\ a,r”). 

The error terms on the right hand-side are written as 

Fpui(90r-, dean = | 5|,a”), 

¥p,,fidga^-\dern = -^r^dga^-^ _ ^ Fp„i(9|,a-i, r"), 

and we then apply the L^-estimates (6.7) and (6.8). This gives contraction of the iteration scheme in 
C([0,Tj;L^(M'^)) and, passing to the limit, we obtain a solution o G L°°([0,Tj;r*'(M'^)) to the Cauchy 
problem. Continuity in r^(M'^) is recovered by the same arguments as in [BGS07, chapter 10], using the 
norm in r^(]R'^) rather than the 77^-norm. We feel free to skip the details. □ 
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We do not claim that the regularity assumption in Theorem 6.4 is optimal as far as local existence of 
smooth solutions is concerned, but it is sufficient for our purpose. The global existence of weak and/or 
smooth solutions is, of course, a wide open problem. Numerical simulations reported in [MR84] tend 
to indicate that finite time breakdown of smooth solutions should be expected, and Proposition 6.2 is 
clearly a first step towards a rigorous justification of this fact. We postpone the study of such finite time 
breakdown to a future work. 

6.4 Construction of the leading profile 

Corollary 6.5 below is based on Theorem 6.4 and is entirely similar to Corollary 3.5 for the wavetrain 
problem. The only difference is that we restrict here to some finite regularity since the estimate provided 
by Theorem 6.4 is not tame. Let us recall that the smoothness assumption for the source term G in (1.1) 
was made in Theorem 1.11. 

Corollary 6.5. With Ko,Ki G N and G as in Theorem 1.11, there exists 0 < T < Tq, and there exists a 
unique 

a G nfioC'((-oo,r];r^i-i-^(M^^,)), 

solution to (6.5), with a\t<co = 0. Furthermore, the integral of a with respect to the variable 0 G M vanishes. 

Proof. The proof is rather straightforward. Since G G C^((— oo,T];r'^i(M'^)), Theorem 6.4 yields a unique 
solution a G C°((—oo,T];r^i“^(M'^)) to (6.5), with o|t<o = 0. Furthermore, (6.5) automatically yields 
a G C^((—oo,T];r^i“^(M'^)) thanks to Proposition 6.2 and the fact that is an algebra. Then 

the standard bootstrap argument yields 

a G nfioC'((-oo,r];r^i-i-^(M^,,)), 

by differentiating (6.5) sufficiently many times with respect to time. 

Let us observe that the requirement Ki — Kq — 1 > 2A:o + 1 in Theorem 1.11 is used here to apply 
Proposition 6.2 for the term i90pQpui[a, a] in (6.1) and its successive time derivatives. 

Thanks to the property a G a is integrable with respect to 0 G M, and integration of (6.11) 

yields 

^Lop tl — 0 ) o|t<0 — 0 ) 

where a denotes the integral of a with respect to 6. Hence a is identically zero. □ 

Corollary 6.6. Up to restricting T > 0 in Corollary 6.5, for all m G X, there exists a unique solution 

am G X M,)) 

to (5.10) with am\t<o = 0 and am\xa=o = where the real number Zm is defined by Cm = ^mi^m- 

Furthermore, each am has zero integral with respect to the variable 9m G K- 

Proof. From Corollary 6.5, we get a G F^°((— oo,T]i xR^ g), with Kq > l + {d+l)/2 (here we use Ki — 1> 
Kq). Then we solve the Burgers equation (5.10) with prescribed boundary condition am\xd=o = a. The 
theory is similar to that in the standard Sobolev space , and we feel free to use well-posedness in this 
weighted Sobolev space framework. This yields a solution 

am G n^l°oC^((-oo,r];r^°-^(R^ X Re)) , 
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to (5.10) that vanishes in the past. Integration of (5.10) with respect to 6^ shows that the integral of cr„ 
with respect to the variable 0m € M satisfies 


^x^m ~ ^ j ^mUd=0 — 0 , 


and therefore vanishes. 


□ 


Remark 6.7. We observe that there has been a rather big loss of regularity in passing from the trace 
a to the interior functions am- This is due to the following fact: the trace a is obtained by solving a 
Cauchy problem, where tangential regularity with respect to the time slices {t = constant} is propagated. 
However, construeting a local in time smooth solution to (5.10) requires a full regularity for the trace 
a, that is, regularity of both {y, 9) and time partial derivatives. This is the reason why we also need to 
control time derivatives of a, which means controlling time derivatives o/50oQpui[o; «]• In Corollary 6.5, 
we have considered sufficiently smooth initial data so that such time derivatives are controlled by an easy 
application of Proposition 6.2. Again, we do not aim at an optimal regularity result. 


6.5 Extension to more general N x N systems 

The extension of the definitions {Vf,Vh, Ep, Eq, Roo, etc.) and resnlts for pnlses in the strictly 
hyperbolic 3x3 case to the strictly hyperbolic N x N case is straightforward. We first show that the 
leading profile Uq reads 

Ko{t,x,9o,^d) = ^ o-m(t,a:,0o + , 

mGX 

with functions dm that are expected to decay at infinity with respect to 9m- Moreover, the (Tm’s satisfy 
(5.10) in the domain {xd > 0}, and (5.4)(a) yields 

Vm € X, am{t,y,0,9o)rm = a{t,y,9o) Cm , 

for a suitable function a. 

The existence of a bounded corrector U 2 solution to (5.2) (c) implies that the first corrector reads 

r +00 

Ui{t,x,9o,^d) = ^ - imEo{t,x,9o + ui^{^d- X)-,X)dX rm 

r- f'~\ 


m£0 


+ V frm(t,a;,0o + WmCrf) - [ Im Eo{t,y,0,9o + UmiU - X), X) dx] rm- 
m&X V J 

Plugging these expressions in (5.4) (b), we derive the following amplitude equation that governs the evo¬ 
lution of a on the boundary; 

(6.11) vde^io?) - ^Lopfl + SeoQpuifajo] =hd0^G, 

with V and Xlop defined in (2.23), and Qpui defined by: 

( 6 . 12 ) 

^ + 00 

Qpul (^ 0 ) *— ^ ^ ^ ^ ^m^mi,m2 j ^(^0 —m) ^(^0 “1“ (^m 2 —m.) 


meO mi<m2 

mi ,m 2 GX 


r +00 

T ^ ^ bB rm ^ Im Em2,mi J ®(0O T (hlmi him) X) 50^(2(00 T (iiim2 him) X) d^ , 


mGO mi<m2 

mi ,m 2 GX 
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with s-s in (2.25). 


As in the wavetrain case some care is needed in the extension to hyperbolic systems with constant 
multiplicities. Instead of introducing bases {im,k}-, m = 1,..., M, k = 1,..., Vkm of ioff nnd right 

eigenvectors, we now define the spaces Vf and Vh and operators Ep, Eq, R^c as follows. 

Definition 6.8. (a) Define the set of “constituent functions” C = where is the set o/M^- 

valued functions of (t, x, 9q, of the form E{t, x, 9 q + where E{t, x, 9m) is and decays with its 

first-order partials at the rate 0{{9)~‘^) uniformly with respect to {t,x). Setting Vc '■= ®m=iC'^j we can 
write any F€Vc as F = Y)m=i with F'^ e . 

(b) Define Vf to be the space of -valued functions of {t,x,9o,f,d) that may be written as a finite 
sum of terms of the form 

F, Tp{K,L,M), 

where F, G,... ,M lie in Vc, Ba ■ x 1^^ —^ is any bilinear map, and Tp : x x —>■ 

is any trilinear map. 

(c) Define a transversal interaction integral to be a function of the form 

r+co 

(6.13) limit,X,9o,u) = - / Ad{t))-^Q^B-,{F\G”^){t,x,9o+ui^{U-s),s)ds, 

hd 

where i,i,m lie in {!,... ,M} and are mutually distinct, By : x —>■ is any bilinear map, and 

F^l G C^, G"* G are required to decay with their first-order partials at the rate 0{{9)~^) uniformly with 
respect to {t,x). 

(d) Define Vh to be the space of -valued functions of {t,x,9Q,f,d) that may be written as the sum 
of an element of Vf plus a finite sum of terms of the form 

(6.14) Il^{t,x,9o,U) or BsiAfi Jl,j{t,x,9o,U)), 

where and are transversal interaction integrals, Bs '. x —)■ is any bilinear map, and 

A^ G is required to decay with its first-order partials at the rate 0{{9)~^) uniformly with respect to 
{t,x). 

(e) For H G Vh we can write H = Hf + Hj, where Hf € Vf and Hj ^Vf is a finite sum of terms of 
the form (6.14). The “constituent functions” of H include those of Hf as well as the functions like Afi 
F^, G™ as in (6.13), (6.14) which constitute Hp 

(f) With these definitions ofVF, Vh, and “constituent functions”, the subspaces V^ and Vp may be 
defined just as in Definition 5.8. 

Definition 6.9 {Ep,Eq,Roo)- For H G Vh, define averaging operators 


M 1 f'^ 

EQH{t,x,9o,^d) lixa - / Qm H{x,9o + uimiCd - s), s) ds , 

T-^oo 1 Jo 


m=l 

M 


Ep H{t,x,9o,f,d) := 


m=l 


6*0 +Wm(Cd - s),s)ds. 


For H E Vh such that Eq H = 0, define the solution operator 


^ f‘-\-00 

RooH{t,x,9o,Cd) ■■=-F / Ad{0)~^ QmH{t,x,9o + Um{f.d-s),s)ds . 

m=l Hi 
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With these dehnitions, Propositions 5.6 and 5.7 hold with the obvious minor changes. For example, 
in Proposition 5.7(c), we now have 


M 

EpU =^W\ where G C™ and = 7/™ . 

m=l 


The construction of prohles is carried out assuming U 2 G and that Uq, U.\ satisfy (5.9), where (5.9) (a) 
is modified to 


M 

m=l 


U^{t,X,er. 


dOr, 


= 0 . 


Again the interior leading profile equations for the QmU^ take the form (4.17), which allows the moment 
zero property of the to be deduced from that of the amplitude a as before. The solvability in T^ 
spaces of (4.17) with boundary conditions (4.20) follows from Lemma 4.2 via estimates proved in the 
standard way.^^ One finds as before that Uq is purely incoming (4.19). 

The consequence (4.16) of the conservative structure assumption is used again in step 2 of the prohle 
construction when performing the integral that now replaces the integral on + 00 [ of the second line 
of (5.14). This integral now reads 



Adior^ 




s), s) ds. 


The replacements for the other lines of (5.14) are treated essentially as before; the last line now yields 
transversal interaction integrals of the form (6.13). 

Parallel to (4.21) the nonlocal operator in the equation for a now has the form 


Qpui[®) o](^0) • 

/ +00 

(90O®(^O + iklmi —m) ®(^0 + (film ,2 —m) 


meO mi<m2 

mi ,m2 


r+oc 

+ E E bP Ad{0) QmPm2,mi I ®(^0 + (l 2 mi ^m) P (—m2 Llm) ^) , 

-*7 0 


mGO mi<m2 

mi ,m2 GX 


with £'mi,m 2 3'S in (2.25). 

Repetition of step 6 of the profile construction yields T/q, Ui, U 2 with the same regularity and decay 
properties as before. 


^^Such an argument is given in detail in Propositions 3.5 and 3.6 of [CW13]. 
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Part III 

Appendices 


A Example: the two-dimensional isentropic Euler equations 

In this first appendix, we discuss how our main results apply to the two-dimensional Euler equations in 
a hxed half-plane. Once again, we refer to [MR83, AM87, WY14] for the derivation of such amplified 
high frequency expansions in various free boundary problems. Our discussion here will mainly deal with 
the verification of the small divisor condition, that is, Assumption 1.9. In quasilinear form, the isentropic 
Euler equations read 


(A.l) 


dtv + ui div + U2 d2V - V {diui + CI 2 U 2 ) = 0, 


dtui -I- ui diui + U 2 d 2 Ui — 
dtU2 -I- Ul dlU2 + U2 d2U2 — 


c{vy 

V ^ 
c{vY 


div = 0. 
d 2 V = 0. 


where u > 0 denotes the specific volume of the fluid, c{v) > 0 denotes the sound speed and (ni,U 2 ) the 
velocity field. We consider a hxed reference volume y_, and a hxed velocity held (0,u) with 


u>0, 0 < u < c := c(v) 


which corresponds to an incoming subsonic huid. The above theory applies when looking for WKB 
solutions to (A.l) of the form 


0 I +£ (t,x, 




n>0 

provided that the linearization of (A.l) (with appropriate boundary conditions) at the reference state 
(v,0,u} satishes all assumptions of Section 1. Let us therefore consider the linearization of (A.l) at 
(u, 0,u), which corresponds, in the notation of Section 1, to 

0 -V 0\ 

Ai(0) := I -cVu 0 0 , A2(0) := 

0 0 0 / 



^U 2 



— V 


1 := I {—2cc'lv + c^lv‘^)v 

Ul 

0 

1 V "" 0 “ 

0 

UiJ 

\ f 

0 

-v\ 

:= 0 

U2 

0 

/ \(—2 cc'/u-|-c^/u^) u 

0 

U 2 ) 


where c' stands for c'(v). The operator dt + Ai(0) di + ^ 2 ( 0 ) 82 is strictly hyperbolic with eigenvalues 
MUi,^2) '■= U^2 - , ■\2Ul,^2) '■= U^2 , -^ 3 (^ 1 , ^ 2 ) '■= 11^2 + C\l + ?! , 
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so Assumptions 1.1 and 1.2 are satisfied. There are two incoming characteristics and one outgoing 
characteristic, so the matrix B encoding the boundary conditions for the linearized equations should be 
a 2 X 3 matrix of maximal rank. One possible choice for B is made precise below. The hyperbolic region 
B is given by 

H = |(r, ?7) £ M X M/ |r| > a/c^ — u? |r?|| . 

We focus on the verification of Assumption 1.9. For concreteness, let (r, r/) £ B with r > 0 and ry > 0, 
and consider the boundary phase 

ipo{t,Xl) :=Tt + T]Xl. 

The three (distinct) eigenvalues of A{t, ry) are 


UT — — (c^ — U^)rj^ UT + Cy^r^ — (c^ — U^) T]^ 

• O O ? ^2 • 






Ws :=-, 

u 


and they satisfy 

r + Ai(r 7 ,a;i) = r + Xi{r],uj 2 ) = r + X 2 {r],L 03 ) = 0. 
The associated group velocities are 


1 


vi := 


2 

-C T] 


1 


+ UUJI — (c^ — ]B) rfj ’ ' t + UOJ2 \—c — (c^ — v'^ 


2 

-c rj 


, V3 := 


hence the phase (/?2 is outgoing while (pi,^P 3 are incoming (as in the framework of Paragraph 2.2). The 
nonresonance condition in Assumption 1.9, that is. 


detL(d(a-$))/O, 


holds if and only if the (dimensionless) quantity 

is not a rational number, as follows from a straightforward calculation. We thus introduce a positive 
irrational parameter C defined by 

(A.2) — a/t^ — (c^ — rf = . 

Our choice of parameters gives C £ (0, 1). For all aijOs £ Z \ {0}, we compute 

det L{d{ai ipi + as y^s)) = -ai as (r + uuji) [2 ai (ry^ + coi W 3 ) + as (ry^ + a;|)] 

= ai as (r + u wi) ^ ~ _ (C - 1) (2 ai C + as (C + 1)) • 

It appears from the latter expression that the verification of Assumption 1.9 only depends on the arithmetic 
properties of the parameter in (A. 2). In particular, we have the following result. 

Lemma A.l. For all 7 > 1, there exists a set of full measure M-y in (0,1) such that for all M £ M-y 
and c > 0, if u = M c and r = cry, then Assumption 1.9 is satisfied for some constant c > 0. 
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Of course, one could also fix the parameters u, c and make r vary in the hyperbolic region, and obtain 
a similar result (meaning that Assumption 1.9 would be satisfied except possibly for r in a negligible set). 

Proof. Choosing t = crj with 77 > 0, one has (t,77) € PL, and the parameter in (A.2) equals with 
M := u/c. Using the above expression for the determinant of L(d(ai ifi + 03 993 )), we see that the small 
divisors condition of Assumption 1.9 will be satished provided that there exists a constant c > 0 such 
that 


(A.3) V(p,g) G \ {(0,0)}, \p + qM^\>c\{p,q)\'^. 

For c G (0,1) and 7 > 1, let us define 

A/; := |777 G ( 0 , 1 )/3 (p,^) G , q^O and \p + q m\ < 


The set A/} is the countable union, as p varies in Z and q in Z*, of intervals of width at most 2 c/dgi \{p, <?)|"’'). 
Hence the Lebesgue measure of A/} is 0(c), which means that the intersection nc>oA/’c is negligible. 
Consequently, for every fixed 7 > 1 , the set of parameters such that (A.3) is satished for some 
constant c > 0 has full measure 1. The claim of Lemma A.l follows. □ 


We assume from now on that r ,77 are hxed such that t = cr], and the Mach number M is chosen in 
such a way that Assumption 1.9 is satished (which is some kind of a generic condition on M). Then we 
compute 

/1 + M 2 \ 


and 



r2 := 


1 -M 2 - 
c 

2u 

V 1 - M2 J 



h := 


1 


1 


2vu 2uc 


0 


^2 := - 


1 


1 + M2 V 2vu 


1 + m2 1 - M^ 1 


2uc 


1 


0 


1 1 


uc 


1 + A /2 

We now make the choice of B precise. As in [CGW14, Appendix B], we choose 


B := 


0 7 ; 0 

u 0 V 


which does not have any physical interpretation but makes the following calculations rather easy. The 
reader can check that Assumptions 1.3 and 1.6 are satished, and we can choose e := ri — r 3 as the vector 
that spans ker B n E^(t, 77 ). The one-dimensional space B ]E^(r, rf) can be written as the orthogonal of the 
kernel of the linear form b:= {u, — c). We can then compute the bilinear Fourier multiplier Qper dehned 
in ( 2 . 20 ), and get: 


Qper [®) o] 


= VUC 


E 

fcGZ 


E 

/ciH-/c 3 =fc, 


1 + 


2 ( 1 - m 2 ) ki 

h + {2h + h)M^ 




J2i7r kOo/S 
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whose kernel is unbounded and depends, as expected, on the arithmetic properties of 

To conclude this Appendix, let us observe that in three space dimensions, the isentropic Euler equations 
are no longer strictly hyperbolic but they enjoy a conservative structure (in the physical variables p,pu). 
Similarly, the full Euler equations with the energy conservation law also have a conservative structure. 
Hence Assumption 1.2 is satisfied, and one can perform a similar derivation as above for the leading 
amplitude equation. 


B Formal derivation of the large period limit: from wavetrains to 
pulses 


B.l The large period limit of the amplitude equation (2.19) 

In this Appendix, we study the relationship between the quadratic operators arising in the leading ampli¬ 
tude equations for the wavetrains and pulses problems which we have considered. Eor the sake of clarity, 
we focus on the easiest case N = 3, p = 2, that was considered in paragraph 2.2 and section 5. The 
Eourier multiplier Qper is then defined by (2.20), while the bilinear operator Qpui is defined by (5.20). 
Erom these expressions, we can decompose both operators as 

Qper ^ ] = {bB r 2 ) h Ei^sWpevidea, a] + (6Hr2) ^2-^3,1 IFper[a, dgo] , 

Qpni[a,a] = {bB r 2 ) £2 Ei^^Wpnildga, a] + (6H r2) 4-E' 3,1 ^puifa, , 


with (observe the slightly different normalizations with respect to (3.3) and (6.2)): 


(B.l) 


iQ 


( 


\ 


^Per[^>^](^) ■“ 2 7 r ^ ^ 




kez 

n+oo 


k^+k3=k, (-1 “ - 2 ) + ^3 - hf2) 

\ ^1^37^0 


^2 ITT kO/Q 


^ p+00 

(B.2) Fpui[u, u](0) := / ^t(0 + (w;^ — hl2)'*)'^^(^ + (fk3 “ ^ 12 ) ®) ds • 

Jo 


In (B.l), both functions u, v are assumed to be 0-periodic and u^, Vk stand for their A:-th Eourier coefficient, 
while in (B.2), u, v are defined on M and have sufficiently fast decay at infinity (so that the integral makes 
sense). 

Our goal is to explain how one can compute the (formal) limit of Fper when the period 0 becomes 
infinitely large and to make the link with the expression of Fpui in (B.2). The additional variables {t,y) 
play the role of parameters here, so we focus on Fper,Fpui as operators acting on functions that depend 
on a single variable 6. We pick two functions u,v in the Schwartz class 5(M), and define 


V0gM, ue{9) :=+ nQ), ve{0) := v{0 + n0). 

nG'Z n^7j 


The functions ue,ve are 0-periodic and converge, uniformly on compact sets, towards u,v as Q tends to 
infinity. Moreover, the Poisson summation formula gives the Fourier coefficients of uq,vq in terms of the 
Fourier transform of u, v: 


“e(») = 5 E 


0 


u 


fcez 


2k-K 


^2i-Kk9/& 




0 


kez 


2k'K 


^2i'Kk9/Q 


58 







We compute 


l^per [^0)'^^ 0 ] (^) — 


2 7 r 0 


E 


kGl. 

i 47r^ 


E 


fcl+fc3=fc, 

\ 

E E 


u{2 ki 7 r/ 0 ) v{2 /c 3 tt/Q) 

ki (wi - W2) + ks (wg - U2) 


^2i7r kO/Q 


/ 


u{2 ki tt/Q) v{ 2 k^ 7 r/ 0 ) 


^2 2 TT /C 0/© 


^1 ^ 37^0 


The latter expression suggests that Fper[M 0 , U 0 ]( 0 ) is the approximation by a Riemann sum, of the double 
integral^^ 

^ C f _ ^r])v{^-r]) _ 




V (wi - ^ 2 ) + i^-r]) (^3 - ^ 2 ) 


d ?7 d.^. 


Formally, this means that, as 0 tends to inhnity, Fper[rt 0 ,ue] tends towards a function whose Fourier 
transform is given by 


(B.3) 




1 

'2i 


L [ _ 

Jr V (hii 


u{r]) v{C - rf) 


-^2) + {i-Tl) (^ 3 -^ 2 ) 


drj, 


assuming of course that the latter expression makes any sense. We wish to compare (B.3) with the Fourier 
transform of Fpui[u, u], whose expression is given by the following (rigorous!) result. 

Lemma B.l. Let u,v € <S(M). Then (B.2) defines a function ¥p^^[[u,v] G 5(M) whose Fourier transform 
is given by 


(B.4) 


2 1^3 — Wp I 


hl3 ~ ^2 
W3 -iill 


(, V 


^2 -^1 
^ 3-^1 


2i7r 


p.v. 


u{ri) v{f, - T]) 


V (hli - ^ 2 ) + i^-v) {^3 - ^ 2 ) 


dr/. 


where p.v. stands for the principal value of the integral. 

Proof. That Fpui[rt,r;] belongs to <S(M) is done in two steps. One hrst proves differentiability by applying 
the classical differentiation Theorem for integrals with a parameter. Then the formula 


Fpui[M, v]' = Fpui[rt', u] + Fpui[tt, v'] , 


and a straightforward induction argument yields infinite differentiability. Given an integer N, we can 
apply the Peetre inequality and get 

^ /■+00 

(1 + (9^)^ |Fpui[rt,r;]((9)| < C sup((l + t^)^ |u(t)|) / (1 + ((9 + (wp - ufi) sf)^ \u{e + (wp - 0 / 2 ) s)|ds, 

teR Jo 

so (1 + 0^)-^ Fpui[rx, u] is bounded. The previous formula for Fpui[rt,r;]' shows again by induction that all 
functions (1 + Fpujrt, are bounded. We may thus compute the Fourier transform of Fpui[rr, u]. 

^“^It is not obvious at first sight that this double integral makes any sense, but our goal here is to identify formally the 
large period limit so let us pretend that all manipulations on the integrals and limits are valid. 
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Let us first write 
(B.5) 
with 


Fp,iKu](0) = Fi(0) + F2(0), 


^i(^) := o +(Lii-W2)'S)^^(6' +(his-Li2)'S)ds, 

^ 2 ( 0 ):=- / sgn(s)u( 6 » +(wi-W 2 ) 5)^(0 + (^ 3 -W 2 )s)ds. 

JR 

Here sgn denotes the sign function. The Fourier transform of Fi is computed by applying an elementary 
change of variables and the Fubini Theorem; 

Fi(^) = - [ e“*^^ u (0 + (^3 — W 2 ) s) u (0 + (^3 — W 2 ) s) ds d 0 
2 Jrz 

X W2) - H (wi - CJ2) 

-f 

2 Jr 


■ 


his -hli 


u{X)v{Y) 


dXdY 
ks -hfil ' 


The contribution of Fi(^) therefore gives the term in the first line of (B.4). It remains to compute the 
Fourier transform of F 2 . 

The Fubini Theorem gives 

^ 2(0 = ;^ [ sgn{s) ( [ u{9 + - U 2 )s)v{e + {u^-ui 2 )s)de) ds = ^ [ sgn{s) {Us-kVs)iO ds, 

2 jR VJR / dvr 

with 

Us{9) := u{9 + (wi - W 2 ) s) , Vs{9) := v{9 + (^3 - W 2 ) s). 

We thus get 

F2{0 = ^ [ sgn(s) ds, 

dTT Jr V JR / 


which is an expression of the form 


^ f sgn(s)£'g(s)ds, 
Jr 


for a suitable Schwartz function is a parameter here). We can therefore transform the expression of 

-^ 2(0 by using the Fourier transform of the sign function, which yields 


MC) = P-v. [ - S^{r]) drj. 
2 * vr 7 r r? 


The function is given by 

sdv) = 


ks - hli 


u 


r/ + (ws - W2) + LI2) ^ 


U)q — UJ^ 


UJ-\ — UJo 


and a last change of variable gives, as claimed in (B.4): 

F2(0 =-x^P-v. / — 

2 ^ vr Juvikh 


u{r])v{C - T]) 


-^2) + iC-v) (^ 3 -^ 2 ) 


drj. 


□ 
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In view of the decomposition (B.5), and the previous computations of Fourier transforms, we observe 
that the formal limit (B.3) for the Fourier transform of Fper[ue,U 0 ] does not coincide with the Fourier 
transform of Fpui[tt, v] but rather coincides with the Fourier transform of the function F 2 . In other words, 
we have formally obtained that in the large period limit 0 —)■ + 00 , Fper[ue,U 0 ] tends towards 


(B. 6 ) 


1 

2 


sgn(s) u{9 + 


^ 2 ) s) w((9 + (W 3 - W 2 ) s) ds , 


which is a ’’symmetrized” version of the operator Fp^i in (B.2). In particular, for any function a G 5(M), 
the ©-periodic function Fper [ 96*00 > ^ 0 ] formally converges, as 0 tends to infinity, towards 

- I sgn{s) dga{6 + {ui — U 2 ) s) a{9 + — U 2 ) s) ds . 

2 Jr 

Remark B.2. It might he surprising that the small divisor condition in Assumption 1.9 does not seem 
to play any role in the analysis of pulses. However, there remains some kind of trace of this Assumption 
but it is hidden in the functional framework. More precisely, we have shown that for some large enough 
index v, the bilinear operator Fpui satisfies a bound of the form 

||Fpui(96)U, 96)u)||r‘' < C'ljujlr'^ ll^^llr'', 

where T'' is a weighted Sobolev space. Omitting the variables {t, y) for simplicity, it is also possible to 
prove that for all integer u, there exists a bounded seguence {an)neN in such that 

(Onj Fpul(9gClfi, ^ -|-oo. 

In other words, the space F^ is well-suited for studying boundedness o/Fpui and the standard Sobolev space 
H'^ is not. This is not so surprising because on the Fourier side, the kernel of Fpui is unbounded, which is 
a trace of the small divisors in the wavetrains problem. Unboundedness of the kernel requires introducing 
a principal value in Lemma B.l, which relies on some continuity of the integrand. It is therefore not 
surprising that continuity of Fpui holds in a functional space where Fourier transforms have some extra 
regularity properties (which is another way to formulate polynomial decay in the physical variable). 

As repeatedly claimed, this paragraph does not aim at giving a rigorous justification of the large 
period limit. An open question that is raised by the formal arguments given above is: assuming that the 
small divisor condition in Assumption 1.9 holds, proving that the 0-periodic function Fper[u 0 ,u©] does 
indeed converge (for instance, uniformly on compact sets), as 0 tends to infinity, towards (B. 6 ) for u,v in 
the Schwartz class. It is also likely that the amplitude equation (6.1) is ill-posed in H‘'(M.‘^) for any large 
integer u, though it is well-posed in F^. We refer to [BGCTll] for similar ill-posedness issues on nonlocal 
versions of the Burgers equation. 


B.2 What is the correct amplitude equation for Mach stem formation ? 

It is surprising that the formal limit of the amplitude equation (2.19) as 0 tends to infinity does not 
coincide with (5.19) if the operator Qpui is defined by (5.20). In several problems of nonlinear geometric 
optics, namely 

• Quasilinear hyperbolic Cauchy problems [HMR 86 , JMR95, AR02, GR06], 
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• Uniformly stable quasilinear hyperbolic boundary value problems [Wil99, Wil02, CGWll, CW13], 

• Weakly stable semilinear hyperbolic boundary value problems [CGW14, CW14], 

the limit of the amplitude equation for wavetrains does coincide with the amplitude equation for pulses. 
Of course, nonlinear geometric optics has received a more complete description for these three problems 
than for the one we consider here because in all three above problems, exact solutions on a fixed time 
interval have been shown to exist and to be close to approximate WKB solutions (which justifies the 
relevance of the corresponding evolution equation for the leading order amplitudes). 

In view of all available references in the literature, it is therefore natural to wonder whether the 
amplitude equation (5.19), (5.20), which is the one derived in [MR83], does give an accurate description 
for pulse-like solutions to (1.1). Rephrasing the question, is it possible to construct another family of 
approximate solutions that would be determined by solving an amplitude equation that is obtained as the 
large period limit of (2.19), (2.20) ? The answer is yes, but the price to pay is to forget about boundedness 
of the second corrector U 2 - 

Let us quickly review the analysis of the WKB cascade (5.2), (5.4) for pulses. In Step 1 of Paragraph 
5.2, we have derived the expression (5.12) of the leading order profile Uq by assuming that the first 
corrector Ui is bounded. The amplitudes solve the decoupled Burgers equations (5.10) and satisfy 

the trace relation (5.11). In Step 4 of Paragraph 5.2, we have obtained the expression of the component of 
EplAi on r 2 by assuming that the second corrector IA 2 is bounded. However, this boundedness assumption 
is not necessarily supported by a clear physical interpretation, and it might very well be that IA 2 is not 
uniformly bounded in the stretched variables 

One can indeed construct (inhnitely) many families of approximate solutions to (1.1). Let us choose 
a parameter s G [0,1]. Then, given the leading order profile Uq in (5.12) with functions (Ti,(T 3 that are 
expected to decay sufficiently fast at infinity, one can construct the first corrector Ui as a particular 
solution to (5.2)(b). One possible choice is: 


f +00 


Ups : = 


imfFo{t,x,eo + iJyn i^d - X),X)dX 


- ^ I 

m=l '' 

+ {1-8) j frnXo{t,X,0Q+yi^{^d-X),X)dx'^ Tm, 

and the general solution to (5.2) (b) is of the form: 

_ 3 

Ups + ^ T^{t, X, 00 + ^d) rm ■ 


m=l 


The choice in [MR83] corresponds, as in Paragraph 5.2, to s = 1 and causality is invoked to set the 
function equal to zero (see Equation (4.10) in [MR83]). In Paragraph 5.2, we have explained why 
T 2 = 0 can de deduced from the assumption that there exists a bounded corrector solution to (5.2) (c). 

If we do not assume existence of a bounded corrector to (5.2) (c), there does not seem to be a clear 
option for constructing the outgoing part t|, and therefore, given s G [0,1], the corrector we have at our 
disposal satisfies 


Ui^s{t,y,O,0o,O) =-kri- s 



hXo{t,y,0,9o 


W2 X, X) dX rs 


+ {l-s) [ ^ 2 -To(t,y, 0 , 6 'o-W 2 X,X)dXr 2 -b*r 3 . 

J —00 
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Plugging the latter expression in (5.4) (b), applying the row vector b and differentiating with respect to 
6q, we obtain Equation (5.19) with the following new definition of the operator QpuP 

^+oo 

Qpui[a,a](6'o) := (bB r2) £2 s / de^aiOo + (w^ - W2)X)o(6'o + (wg -U2)X)dX 

Jo 


f + OO 


+ [bB r2) £2 s / fl(^o + (lili ~ ^2) ' 96 )oo(^o + (^3 “ 1^2) 

Jo 

- (6Sr2)£2 ^ 1,3 (1 - s) [ dega{0o + - U 2 ) X)a{9o + - ^ 2 ) 

J —00 

(B.7) -(6Br2)£2^3,i(l-5) f a{eo + {^^-ui 2 )X)dg,a{eo + {ui^-^ 2 )X)dX . 

J —00 

Unlike all other terms in (5.19), the bilinear term Qpui does depend on the value of s that is chosen for 
constructing the particular solution Ui^s to (5.2) (b), and therefore invoking causality to discard the r| 
term has an impact on the leading order amplitude equation (5.19). The choice s = 1/2 is the only that 
is compatible with the large period limit 0 —>■ + 00 . 

Our overall assessment is the following. In the wavetrain problem, the construction of approximate 
WKB solutions in Part I is well-understood. In particular Theorem 1.10 shows that the full WKB cascade 
has a unique solution (7/„)„gN in a- suitable functional space. There remains of course the problem of 
understanding the lifespan of exact solutions to ( 1 . 1 ) and whether approximate and exact solutions are 
close to each other. In the pulse problem, the situation is worse since the construction of approximate 
solutions in Part II may be questionable. There are basically two options (the infinitely many other ones 
seem to be a mathematical artifact). Either one hopes that the expansion of exact solutions will yield 
a second corrector U2 that is bounded, and in that case the appropriate leading amplitude equation is 
given, as in [MR83], by (5.19), (5.20). Or one rather expects that the amplitude equation for pulses 
should coincide with the large period limit of the corresponding wavetrain equation, and in that case the 
appropriate leading amplitude equation is still given by (5.19) but with the bilinear operator Qpui as in 
(B.7) with s = i-j’l. The two corresponding leading profiles differ in each of the two options. Without any 
precise understanding of the behavior of exact solutions, deciding between the two possible expansions 
seems hardly possible. 


C Some remarks on the resonant case 

In this Appendix, we explain why the analysis in Part I breaks down when resonances occur. Let us 
consider for simplicity the 3x3 strictly hyperbolic case of Paragraph 2.2. A resonance corresponds to 
the existence of a triplet (ni,n 2 ,n 3 ) G Z^, with ni re 2 n 3 / 0 , gcd (ni,n 2 ,n 3 ) = 1 , and 


(C.l) 


m ipx = n2^2 + ns ipz ■ 


Incoming waves associated with the phases may then interact to produce an outgoing wave for the 

phase p2- More precisely, the source term G{t,y,ipo/e) in (1.1) is expected to give rise to incoming 
waves associated with the phases pi, (^3 of amplitude 0(e). Then the products Aj{us) djU^ in (1.1) ignites 
oscillations associated with the outgoing phase p 2 of amplitude^^ 0(e). These y 52 -oscillations are then 


^^This is a major scaling difference with our previous work [CGW14] where nonlinear interaction was due to a zero order 
term so the outgoing oscillations produced by resonances had amplitude 0{e^), instead of Ole) here. 
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expected to be amplified when reflected on the boundary, as in the linear analysis of [CGIO], but this 
would give rise to , yjs-oscillations of amplitude 0(1) and would thus completely ruin the ansatz (2.5). 
Resonances are therefore expected to produce dramatically different phenomena from the ones that are 
studied here. 

Another hint that resonances should make the ansatz (2.5) break down is to try to solve the WKB 
cascade (4.1), (4.3) when the resonance (C.l) occurs between two incoming phases (^i, ^33 and one outgoing 
phase ip 2 - Let us recall that in the 3x3 strictly hyperbolic case considered in Paragraph 2.2, B is 
a 2 X 3 matrix of maximal rank. Its kernel has dimension 1 and is therefore spanned by the vector 
e = ei + 63 G ]E^(r, 77 ). In other words, we can choose a vector e € E*(t, r/) such that 

= Span (B e, B r 2 ). 


We now follow the analysis of Paragraph 2.2 and try to analyze the WKB cascade in the presence of 
a resonance. Equation (2.6) (a) shows again that the leading prohle can be decomposed as 

3 

l(o(t,X,0i,02,93) =llo(t,x) + ^ am{t,X,0m)rm- 

m=l 

Then Equation (2.6)(b) shows that the mean value Hq satishes (2.10), (2.11), and therefore vanishes. The 
interior equation satisfied by each am exhibits the resonance between the phases, see [Raul2]. Namely, 
the fJm’s satisfy the coupled Burgers-type equations 


(C.2) 


dtcri + Vi • VxCTi + Cl ai de^^ai = Bi{a 2 ,a 3 ), 

< dta 2 + V2 • Vx(T 2 + C2 a 2 502 fJ2 = B 2 {ai,a 3 ) , 

dta3 + V3 • Va;CT3 -b C3 0-3 503(73 = B 3 {ai,a 2 ), 


where the constants Cm are defined as in (2.12) and, for instance: 


Bi{a2,a3){t,x,0i) 


2i7r£idj(p3 (dAj(O) - 72)73 

0.^1 ri 


^ Ckn2{(^2){t,x) (A:n3)cfcn3(o'3)(La:)e^*’"^"^®^/® 

fcez* 


+ 


2 f vr djip2 (dAj(O) • 73) r2 
Q£i ri 


^ (A:n2)cfc„2(o'2)(La:)cfc„3((T3)(t,x)e^*’"^"i. 
fcez* 


The definitions of B 2 {ai,a 3 ) and B 3 {ai,a 2 ) are similar, and correspond to the so-called interaction inte¬ 
grals in [CGWll]. 

The boundary conditions for the (T^’s are given by (2.8)(a). Using the basis (e, e) of Span (ri, rs), we 
decompose 

^o(L y, 0, 00 , 00 , 0o) = a{t, y, 0o) e -b a{t, y, 0o) e -b a 2 {t, y, 0, 0o) r 2 , 

so (2.8) (a) gives 

0 = 0, a2\xa=0 = 0 . 

However, this boundary condition on a 2 does not seem to be compatible with (C.2) because 02 satisfies 
an outgoing transport equation with a nonzero source term (it is proved in [CGWll] that the operators 
Bi,B 2 ,B 3 in (C.2) act like semilinear terms and do not contribute to the leading order part of the 
differential operators in (C.2)). Even if we manage to isolate an amplitude equation for determining the 
traces of the incoming amplitudes ( 71 ,( 73 , the overall cascade seems to give rise to an overdetermined 
problem for Uo- 
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We do not claim of course that the above arguments give a rigorous justification of the ill-posedness 
of (4.1), (4.3) when a resonance occurs, but it clearly indicates that the ansatz (2.5) does not seem to be 
appropriate any longer. 

When there is only one incoming phase, all the above discussion becomes irrelevant, and this may 
again be one reason why any discussion on resonances is absent from [AMS7] or [WY13, WY14]. 
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